Defense  Nuclear  Agency 
Alexandria,  VA  22310-3398 


DNA-TR-96-19 


Load-Damage  Relationships  for  Chemical 
Submunitions 


James  D.  Colton,  et  al. 

SRI  International 

333  Ravenswood  Avenue 

Menlo  Park,  CA  94025-3434 


October  1996 


Technical  Report 


CONTRACT  No.  DNA  001 -93-C-01 04 

a 


Approved  for  public  release; 
distribution  is  unlimited. 


19961022  030 


DnC  QUALETY  IN8PSCTED  1 


DESTRUCTION  NOTICE: 


Destroy  this  report  when  it  is  no  longer  needed. 

Do  not  return  to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  SPECIAL  WEAPONS 
AGENCY,  ATTN:  CSTI,  6801TELEGRAPH  ROAD. 
ALEXANDRIA,  VA  22310-3398,  IF  YOUR  ADDRESS  IS 
INCORRECT,  IF  YOU  WISH  IT  DELETED  FROM  THE 
DISTRIBUTION  LIST,  OR  IF  THE  ADDRESSEE  IS  NO 
LONGER  EMPLOYED  BY  YOUR  ORGANIZATION. 


CUT  HERE  AND  RETURN 


DISTRIBUTION  LIST  UPDATE 


This  mailer  is  provided  to  enable  DSWA  to  maintain  current  distribution  lists  for  reports.  (We  would 
appreciate  vour  providing  the  requested  information.) 

□  Add  the  individual  listed  to  your  distribution  list. 

□  Delete  the  cited  organization/individual. 

□  Change  of  address. 

NAME:  - - 


NOTE: 

Please  return  the  mailing  label  from  the 
document  so  that  any  addHIons,  changes, 
corrections  or  deletions  can  be  made  easily. 
For  distribution  cancellation  or  more 
information  call  DSWA/IMAS  (703)  325-1036. 


ORGANIZATION:  - 

OLD  ADDRESS 


CURRENT  ADDRESS 


TELEPHONE  NUMBER:  i _ I 


DSWA  PUBLICATION  NUMBER/TITLE 


CHANGES/DELETIONS/ADDITIONS,  etc.) 

{Attach  Sheet  it  more  Space  is  Required) 


DSWA  OR  OTHER  GOVERNMENT  CONTRACT  NUMBER: 


CERTIFICATION  OF  NEED-TO-KNOW  BY  GOVERNMENT  SPONSOR  (if  other  than  DSWA): 

SPONSORING  ORGANIZATION:  - - - - 

CONTRACTING  OFFICER  OR  REPRESENTATIVE:  - - 


SIGNATURE: 


DEFENSE  SPECIAL  WEAPONS  AGENCY 
ATTN:  IMAS 

6801  TELEGRAPH  ROAD 
ALEXANDRIA,  VA  22310-3398 


DEFENSE  SPECIAL  WEAPONS  AGENCY 
ATTN:  IMAS 

6801  TELEGRAPH  ROAD 
ALEXANDRIA.  VA  22310-3398 


Form  Approved 
0MB  No.  0704-0188 


REPORT  DOCUMENTATION  PAGE 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response  including  the  time  for  reviewing  instnjctions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services  Directorate  for  information  Operations  and  Reports.  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

961001 


4.  TITLE  AND  SUBTITLE 

Load-Damage  Relationships  for  chemical  Submunitions 


6.  AUTHOR(S) 

James  D.  Colton,  Curtis  M.  Romander,  Steven  W.  Kirkpatrick,  Alexander  L,  Florence, 
and  Jacques  H.  Giovanola 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025-3434 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
Defense  Special  Weapons  Agency 
6801  Telegraph  Road 
Alexandria,  VA  22310-3398 
WELR/Hunter 


3.  REPORT  TYPE  AND  DATES  COVERED 
Technical  930520-950131 


5.  FUNDING  NUMBERS 
C  -DNA001-93-C-0104 
PE  -63216C,63216C,63214C 
PR  -  SB,  AJ,  AJ 
TA  -  AE,  AE,  AF 
WU  -  DH338260 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

PYU-4727 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

DNA-TR-96-19 


11.  SUPPLEMENTARY  NOTES 

This  work  was  sponsored  by  the  Ballistic  Missile  Defense  Office  (BMDO)  and  managed  and  executed  by  the  Defense  Special 
Weapons  Agency  under  RDT&E  RMC  Codes  B7668D  SB  AE  00002  7010A  AJ  25904D,  B7668D  AJ  AE  00002  7010A 
25904D,  and  B7666D  AJ  AF  00012  7010A  25904D. 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


1 3.  ABSTRACT  {Maximum  200  words) 

We  determined  the  load-damage  relationship  for  individual  chemical  submunitions  to  link  (1)  the  loads  on  individual 
submunitions  in  a  complete  target  that  can  be  measured  or  calculated  and  (2)  the  damage  to  individual  submunitions  that  is 
more  difficult  to  determine.  Impact  experiments  and  calculations  were  conducted  on  single  Chemical  Submunitions  1  and  2 
(nominally  95%  full  of  water)  at  full-,  half-,  and  quarter-scale,  at  velocities  up  to  330  m/s  to  simulate  impacts  away  from  the 
direct  impact  area  in  a  complete  target  impacted  by  a  hit-to-kill  vehicle.  The  damage  mechanisms  were  similar  to  those  pro¬ 
duced  in  sled  tests  of  complete  targets:  failure  of  the  burst  diaphragm,  removal  of  the  fuze  plug,  fracture  where  the  case  is 
welded  to  the  top  and  axial  fractures  in  the  casing.  The  load-damage  relationship  was  quantified  in  terms  of  pressure-impulse  or 
PI  curves.  We  also  determined  that  submunitions  could  be  ejected  at  velocities  of  about  65  m/s  with  no  damage  to  prevent  their 
functioning. 
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fuze  removal,  are  the  same  at  botli  scales.  Failure  of  the  weld  did  not  occur  in  the  quarter-scale  because  of  intentional  differ¬ 
ences  in  construction.  Fractures  in  the  cylindrical  casing  of  the  submunitions,  less  frequently  observed  than  diaphragm  or  fiize 
failure,  were  more  ductile  and  hence  required  greater  equivalent  loads  to  produce  in  quarter-scale  than  in  full-scale. 
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EXECUTIVE  SUMMARY 


Missiles  with  payloads  comprising  arrays  of  chemical  submunitions  are  a  key  target  in  the  Theater 
Missile  Defense  Program.  To  determine  the  respcmse  of  such  payloads  to  hit-to-kill  vehicles,  other 
organizations  have  conducted  sled  tests  at  impact  velocites  up  to  2  knVs,  to  observe  submunition  response, 
and  have  perfOTined  computer  code  calculations  of  the  response.  However,  neith^  approach  allows  a 
detailed  investigation  into  the  load  and  response  of  individual  submunitions.  Submuntions  that  are  away 
from  the  impact  area,  where  the  secondary  impact  velocities  fiom  fragments  or  neighboring  submunitions 
are  hundreds  of  meters  per  second,  are  of  particular  interest  because  the  damage  to  those  submunittMis 
strongly  affects  the  damage  to  the  overall  array.  A  determination  of  the  load-damage  relationships  for 
chemical  submunitions  can  link  measiuements  of  response  made  in  sled  tests  and  calculated  in  the  most 
sophisticated  codes  to  final  damage  to  individual  submunitions.  Load-damage  relationships  are  also 
important  in  performing  lethality  assessments  where  damage,  which  must  be  assessed  for  a  range  of 
parameters,  cannot  be  determined  in  sled  tests  or  with  computer  code  calculations. 

The  primary  objective  of  the  research  reported  here  was  to  determine  the  load-damage  relationships 
for  chemical  submunitions.  Secondary  objectives  were  to  determine  the  maximum  velocity  with  which  a 
surviving  submunition  can  be  ejected  from  the  submunition  array  (for  post-engagement  analysis),  to 
provide  data  for  computer  code  validation  of  submunition  response,  and  to  evaluate  scaling  of  chemical 
submunitions  to  impact  response. 

Our  approach  consisted  of  the  following  steps: 

•  Identify  the  damage  mechanisms  produced  in  the  sled  tests 

•  Perform  experiments  and  analyses  to  investigate  the  damage  mechanisms  in  detail 

•  Based  on  results  of  the  above,  develop  load-damage  relationships  in  terms  of  pressure- 
impulse  relationships 

•  Compare  damage  mechanisms  in  full-  and  quarter-scale  hit-to  kill  tests  on  submunition 
arrays  and  perform  well-controlled  impact  experiments  and  material  property 
measurements  at  different  scales. 

We  conducted  impact  experiments  on  single,  full-scale  Chemical  Submunitions  1  and  2,  and  on 
half-  and  quarter-scale  Chemical  Submunitions  2  at  velocities  up  to  330  m/s.  Most  submunitions  had  their 
volume  filled  to  95%  capacity  with  water  to  simulate  the  chemical  agent  To  study  basic  response 
mechanisms  of  submunitions,  most  tests  were  performed  with  a  gas  gun  to  accelerate  steel  plates  or  actual 
submunition  nKxlels  into  single  submunition  targets.  In  the  impact  experiments,  we  measured  the  impact 
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velocity  and  firequently  used  high-speed  movies  and  high-speed  video  to  observe  the  impacts  and  to 
measure  the  velocity  imparted  to  the  target  We  performed  computer  code  calculations  on  individual 
submuntions  to  plan  and  interpret  the  experiments  and  to  broaden  our  understanding  of  the  response 
phenomena. 

We  conducted  a  series  of  ejqreriments  in  which  single,  full-scale  submuntions  were  impacted  with 
a  plate  to  deteimine  the  locaticm,  mentation,  and  magnitude  of  the  loads  required  to  produce  damage 
similar  to  that  observed  in  the  sled  tests.  Figure  S-1  summarizes  the  results  for  Chemical  Submutution  1 
(similar  experiments  woe  performed  for  Chemical  Submuiution  2).  Impact  of  the  fuze  end  of  the  models 
demonstrated  that  the  fuze  acts  as  a  good  energy  absorber  that  protects  the  body  of  the  submunition. 
Therefore,  nose-on  loads  are  the  largest  loads  needed  to  produce  damage.  The  smallest  loads  to  produce 
damage  are  those  that  occur  on  the  shoulder,  causing  the  weld  that  coruiects  the  top  and  the  main  body  of 
the  submuntion  to  fail.  For  side  impact,  the  damage  mechanism  was  frequently  mpture  of  the  burst 
diaphragm  and/or  removal  of  the  fuze. 

We  performed  plate  impact  experiments  and  computer  code  calculations  to  investigate  the  effect  of 
ullage  amount  and  location.  We  found  that,  for  the  range  of  dynamic  loadings  investigated,  the  failure  of 
the  diaphragm  is  largely  a  dynamic  phenomenon  in  the  fluid.  Therefore,  the  amount  of  ullage  did  not 
significantly  affect  diaphragm  failure.  We  also  found  that  the  location  of  the  ullage  was  not  important  in 
causing  damage  to  the  diaphragm;  the  same  damage  occurred  even  when  the  ullage  was  located  at  the 
diaphragm. 

We  performed  a  series  of  experiments  and  calculations  to  quantify  the  load-damage  relatitmship  in 
terms  of  the  applied  peak  pressure  and  impulse  (PI)  curves.  Figure  S-2  summarizes  the  results  for 
external  loading  that  spans  the  entire  side  of  Submunition  2.  The  figure  illustrates  how  different  types  of 
experiments  can  be  used  to  span  the  load  durations  of  interest  The  figure  also  shows  the  similarity 
between  the  experimental  results  and  the  calculations,  and  how  both  can  be  used  to  construct  PI  curves. 
Techniques  based  on  the  response  characteristics  of  the  diaphragm  were  developed  to  assign  pressures  and 
impulses  to  the  mxjlti-peak,  finite-rise,  waveforms  produced  in  the  fluid.  The  figure  shows  calculated 
critical  internal  pressures  and  impulses  that  correspond  to  the  critical  external  loads.  Similar  data  and 
calculations  provided  additional  points  on  PI  plots  to  characterize  external  spot  loading  (i.e.,  plate  impact) 
on  the  side  and  nose-on  loading  for  Chemical  Submunition  2,  and  side,  spot,  and  nose-on  loading  for 
Chemical  Submunition  1.  The  PI  curves  for  either  external  or  internal  loading  can  be  used  to  link 
measurements  and  calculations  of  loads  produced  in  hit-to-kill  scenarios  with  damage  to  submurutions. 

Based  on  the  mechanisms  observed  in  the  experiments  and  the  code  calculations,  we  developed 
simple  analyses  to  predict  PI  curves  for  diaphragm  rupture  and  fuze  failure.  The  PI  curves  derived  from 
the  simple  analyses  can  be  used  to  predict  response  for  other  parameters  of  interest  For  example,  if  the 
thickness  of  the  submuiution  is  decreased,  the  simple  analyses  can  be  used  to  estimate  PI  curves  for  the 
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300  or  Greater  @  22.5' 


Figure  S-1 .  Impact  velocities  to  produce  failure  in  Chemical  Submunition  1  for  different 
types  of  impacts  (1.9-cm-thick,  8.9-cm-diameter  impactors). 
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Figure  S-2.  PI  curve  for  external  side  load  Chemical  Submuition  2. 
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modified  submunitions.  These  analyses  are  particularly  useful  for  lethality  assessment,  because  a  large 
number  of  loading  c<mditi(Mis  can  be  quickly  evaluated. 

Submuniticm  to  submunition  impacts  were  also  conducted  to  determine  the  maximum  residual 
velocity  that  could  be  imparted  to  a  submunition  that  survived  the  intact  Frcwn  tests  with  die  impacting 
submunition  and  the  impacted  submunition  aligned,  and  tests  with  the  two  submuntions  at  right  angles  to 
each  other,  we  found  that  the  maximum  velocity  that  could  be  imparted  to  a  surviving  submunition  was 
about  65  m/s. 

We  also  investigated  scaling  of  submunition  response.  Available  data  from  full-scale  hit-to-kill 
sled  tests  and  quarter-scale  gun  tests  indicated  that  fragmentation  of  the  submuntitMis  in  the  vicinity  of  the 
impact  appears  to  scale.  Differences  in  the  damage  to  the  submunition  casing  resulted  from  differences  in 
construction  (welds  versus  no  welds — ^not  a  true  scaling  difference)  and  from  diffi^nces  in  fracture 
behavior  between  the  full-  and  the  quarter-scale  material.  The  (scaled)  loads  to  produce  fracture  of  the 
casing  were  greater  at  quarter-scale  than  at  full-scale.  Thus,  if  fracture  of  the  submunition  casing  is  a 
significant  damage  mode  of  interest,  care  must  be  taken  in  interpreting  the  results  of  small  model 
experiments.  When  fragmentation  ot  fracture  of  the  casing  did  not  occur  (which  was  the  case  for  a 
majority  of  the  submuntions  in  the  anay)  and  when  the  burst  diaphragm  static  strength  is  properly  scaled, 
the  damage  mechanisms  and  critical  loads  (notably  diaphragm  and  fuze  failure)  were  similar  in  the  full-  and 
quarter-scale  submunitions. 

Although  the  half-scale  models  had  a  greater  hydrostatic  burst  pressure  to  fail  the  di{q)hragm  than 
the  other  models,  the  failure  mechanisms  were  the  same  as  for  full-scale  and  we  expect  that  properly 
designed  and  construaed  models  would  have  critical  loads  to  cause  failure  very  close  to  that  of  the  full- 
scale. 


We  recommend  that 

•  Additional  experiments  be  carried  out,  particularly  on  Submunition  1,  to  awnplete  the 
critical  load  curves  for  external  and  internal  loads  and  fca*  multiple  impacts. 

•  Simplified  models  be  developed  for  loading  over  a  portion  of  the  side  of  the 
submunition,  for  axial  impact,  and  for  internal  loads  for  loading  other  than  side  area 
loading. 

•  TTie  critical  load  curves  developed  be  integrated  into  lethality  assessment  codes  such  as 
KAPPII. 

•  Additional  experiments  be  conducted  to  refine  the  critical  survival  velocity  imparted  to  an 
impacted  submunition  by  glancing  blows  and  multiple  impacts.  Post-engagement 
an^yses  of  the  results  should  also  be  conducted. 


vn 


PREFACE 


The  work  discussed  in  this  technical  report  was  performed  from  May  1993  through  November 
1994  under  Contract  No.  DNA  001-93-C-0104.  Commander  K.  Hunter  was  the  DNA  COTR.  This 
report  deals  primarily  with  structural  damage  to  chemical  submunitions.  A  companion  technical  report  is 
being  issued  concurrently  entitlted  "Improved  Fragmentation  Algorithms  for  Debris  Environments." 


CONVERSION  TABLE 


£njL?ffi//^ctors  for  U-  S-  Customary  to  pietri(^(gl)  units  of  measurement 
TO  GET  -  —  BY  ^ 


angstrom  LOOOOOO  XE-iu  inet 

atmosphere  (normal)  1.013  25  XE+2  kilo 

bar  1.000000  XE+2  kilo 

bam  1.000  000  X  E  -28  met 

British  thermal  unit  (thermochemical)  1.054  350  XE+3  joul 

calorie  (thermochemical)  4.184000  joul 

cal  (thermochemical)/cm^  4.184000  XE-2  mef 

curie  3.700  000  XE+1  *gil 

degree  (angle)  1.745  329  XE-2  radi 

degree  Fahrenheit  Tk  =  (T®F+459.67)/1.8  deg 

electron  volt  1.60219  XE-19  joul 

erg  1.000  000  X  E  -7  joul 

erg/second  1 *000  000  X  E  -7  wat 

foot  3.048  000  XE-1  met 

foot-pound-force  1.355  818  joul 

gallon  (U.S.  liquid)  3.785  412  X  E  -3  met 

inch  2.540000  XE-2  met 

jerk  1.000  000  X  E  +9  joul 

joule/kilogram  (J/kg)  (radiation  dose  1.0(X)  000  Gra 

absorbed) 

kilotons  4.183  tera 

kip  (1000  Ibf)  4.448  222  XE+3  nev 

kip/inch2(ksi)  6.894757  XE+3  kilc 

ktap 

1.000000  XE+2  (N- 

micron  1 .000  000  X  E  -6  met 

mil  2.540  000  XE-5  met 

mile  (international)  1.609  344  XE+3  met 

ounce  2.834952  XE-2  kilc 

pound-force  (lbs  avoirdupois)  4.448  222  nev 

pound-force  inch  1.129  848  XE-1  nev 

pound-forceAnch  1.751268  XE+2  nev 

pound-force/foot2  4.788  026  XE-2  Idle 

pound-force/inch2  (psi)  6.894757  kilc 

pound-mass  (Ibm  avoirdupois)  4.535  924  XE-1  kilc 

pound-mass-fbot2  (moment  of  inertia)  kilc 

4.214  011  XE-2  (kg 

pound-mass-foot^  kilc 

1.601  846  XE+1  (kg 

rad  (radiation  dose  absorbed)  1.000  000  XE-2  **( 

roentgen  coi 

2.579  760  X  E  -4  (Q 

shake  1 *000  000  X  E  -8  sec 

slug  1.459  390  XE+1  kilc 

torr (mm Hg. 0°  C) _ _ I  1.333  22  XE-1 _ [kilc 

♦The  becquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  =  1  event/s. 

♦♦The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 


-►TO  GET 
DIVIDE 


1.000000 

XE-10 

meters  (m) 

1.013  25 

XE+2 

kilo  pascal  (kPa) 

1.000000 

XE+2 

kilo  pascal  (kPa) 

1.000000 

XE-28 

meter^  (m^) 

1.054  350 

XE+3 

joule  (J) 

4.184000 

joule  (J) 

4.184000 

XE-2 

mega  joule/m^  (MJ/m^) 

3.700  000 

XE+1 

♦giga  becquerel  (GBq) 

1.745  329 

XE-2 

r^an  (rad) 

Tk  =  (T®F+ 459.67)71.8 

degree  kelvin  (K) 

1.602  19 

XE-19 

joule  (J) 

1.000  000 

XE-7 

joule  (J) 

1.000000 

XE-7 

watt  (W) 

3.048  000 

XE-1 

meter  (m) 

1.355  818 

joule  (J) 
meter-*  (m-*) 

3.785  412 

XE-3 

2.540000 

XE-2 

meter  (m) 

1.000  000 

XE+9 

joule  (J) 

1.000  000 

Gray  (Gy) 

4.183 

terajoules 

4.448  222 

XE+3 

newton  (N) 

6.894757 

XE+3 

kilo  pascal  (kPa) 
newton-second/m^ 

1.000000 

XE+2 

(N-s/m^) 

1.000000 

XE-6 

meter  (m) 

2.540  000 

XE-5 

meter  (m) 
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SECTION  1 


INTRODUCTION 


1.1  BACKGROUND. 

Missiles  with  payloads  comprising  arrays  of  chemical  submunitions  are  a  key  target  in  the  Theater 
Missile  Defense  Program.  To  determine  the  response  of  such  payloads  to  hit-to-kill  vehicles,  sled  tests 
have  been  conducted  at  Holloman  Air  Force  Base  as  part  of  a  program  conducted  by  the  U.S.  Army 
Strategic  Defense  Ccmimand  (Mullins  et  al.,  1992).  Sled  tests  have  been  conducted  at  impact  velocites  up 
to  2  km/s  to  observe  submunition  response.  State-of-the-art  computer  codes  have  also  been  used  to 
calculate  the  response  of  the  ccnnplete  submunition  array  (Smith  and  Durham,  1993).  However,  neitho* 
approach  allows  a  detailed  investigation  into  the  load  and  response  of  individual  submunitions,  particularly 
those  that  are  away  from  the  impact  area  where  the  impact  velocities  are  reduced  to  hundreds  of  meters  per 
second.  Such  submunitions  are  large  in  number  and  their  fate  often  determines  the  damage  to  the  overall 
array.  Details  of  the  loads  applied  to  these  submunitions  are  not  available  from  sled  tests,  because  such 
measurements  are  difficult  and  expensive  to  make.  Similarly,  in  calculations  performed  to  date,  external 
loads  on  submuniticais  are  difficult  to  predict  because  of  the  complex  loading  and  response  mechanisms 
for  high  velocity  impacts  and  the  fine  zoning  required  with  current  limitations  on  computing  power. 

Understanding  the  relationship  between  loads  and  damage  for  these  submunitions  is  useful  in 
several  ways.  Load-damage  relationships  are  useful  for  comparing  loads  required  to  produce  damage 
under  different  impact  conditions  and  for  different  submuntions  (i.e.,  to  provide  a  first  order  assessment 
of  target  hardness).  Load-damage  relationships  are  particularly  useful  as  a  framework  to  tie  together 
results  ffcrni  different  tests  and  analyses.  Load-damage  relationships  that  are  simple  in  form  can  be  used  as 
essential  input  into  lethality  assessments  where  damage,  which  must  be  assessed  for  a  range  of 
parameters,  is  not  feasible  to  be  determined  in  sled  tests  or  with  detailed  computer  code  calculations. 
Finally,  a  determination  of  load-damage  relationships  can  link  measurements  made  in  sled  tests  and 
calculated  in  code  calculations  final  damage  to  individual  submunitions.  In  many  sled  tests, 
measurements  have  been  made  of  the  internal  pressure  at  the  center  of  one  or  more  submunitions.  In 
calculations  performed  of  the  response  of  complete  arrays,  it  is  believed  that  the  internal  pressure  can  be 
calculated  mote  accurately  than  many  other  response  parameters.  Thus,  a  determination  of  the  load- 
damage  relaticMiships  for  chemical  submunitions  that  includes  internal  pressure  can  link  final  damage  to 
individual  submunitions  with  measurements  and  calculations  in  complete  submuntions  arrays.  Then,  load- 
damage  relationships  can  be  used  in  performing  lethality  assessments  for  many  scenarios  of  interest. 
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1.2  OBJECTIVES. 

The  primary  objective  of  the  research  reported  here  was  to  determine  the  load-damage  relationships 
for  individual  Chemical  Submunitions  1  and  2  (Figure  1-1)  under  a  variety  of  load  orientations  and 
durations.  The  load-damage  relationships  should  also  relate  externally  applied  loads  to  internal  pressure  to 
the  extent  possible.  Sectxidary  objectives  were  to  determine  the  effects  of  ullage  on  submunition  impact 
response,  to  determine  the  maximum  velocity  with  which  a  surviving  submunition  can  be  ejeaed  from  the 
submunition  array  (for  post-engagement  analysis),  to  provide  data  for  validating  computer  code 
calculations,  and  to  evaluate  scaling  of  chemical  submunitions  to  impact  response. 

1.3  APPROACH. 

To  meet  these  objectives,  we  tt)ok  a  combined  experimental  and  analytical  approach  that  consisted 
of  the  following  steps.  First,  we  examined  the  submunitions  damaged  in  some  of  the  sled  tests  to 
determine  the  dominant  damage  mechanisms.  Second,  we  performed  e^qreriments  and  analyses  to 
investigate  these  mechamsms  in  more  detail,  to  quantify  the  relationship  between  load  and  damage  for 
particular  cases  of  interest,  and  to  determine  the  effect  of  ullage  locations.  We  performed  submunition  to 
submunition  impact  tests  to  determine  the  maximum  velocity  at  which  a  surviving  submunition  can  be 
ejected  from  the  submunition  array.  Next,  we  combined  the  load-damage  relationships  already  determined 
for  particular  cases  and  performed  additional  experiments  and  analysis  so  that  we  could  express  the  load- 
damage  relations  in  terms  of  pressure-impulse  (PI)  characterizations  for  a  variety  of  loading  conditions  of 
rnteresL  Based  on  the  observed  damage  mechanisms,  we  performed  simplified  analyses  that  allow 

prediction  of  the  critical  loads  for  target  parameters  other  than  those  tested  (e.g.,  for  different  submunition 
wall  thicknesses). 

To  investigate  scaling  issues,  we  first  compared  the  damage  observed  in  the  sled  tests  on 
submunition  arrays  conducted  at  Holloman  Air  Force  Base  with  corresponding  damage  observed  in 
quarter-scale  array  tests  conducted  at  the  Arnold  Engineering  Development  Center.  Then,  we  compared 
the  responses  of  tensile  and  fracture  tests  on  material  specimens  from  the  full-,  half,  and  quarter-scale 
submunition  and  of  carefully  controlled  impact,  tests  on  single  fuU-scale  and  corresponding  half-  and 
quarter-scale  submunitions. 
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SECTION  2 

CHEMICAL  SUBMUNITION  IMPACT  RESPONSE 


To  investigate  the  impact  response  of  chemical  submunitions,  we  first  examined  the  damage  modes 
produced  in  the  sled  tests.  We  then  performed  both  experiments  and  analyses  to  investigate  and  quantify 
the  damage  produced  in  these  modes. 

2.1  DAMAGE  MODES  OBSERVED  IN  SLED  TESTS. 

We  visited  Teledyne  Brown  Engineering  in  Huntsville,  where  we  examined  the  damage  to  the 
submunitions  from  the  sled  tests  conducted  at  Holloman  Air  Force  Base.  In  these  tests,  a  mock-up 
simulating  a  hit-to-kill  projectile  was  accelerated  by  a  rocket  sled  and  impacted  a  target  containing 
submunition  threat  simulators  (hereafter  referred  to  as  submunitions).  The  objective  of  our  examination 
was  to  identify  the  damage  modes  in  the  submunitions,  particularly  those  modes  that  would  result  in 
release  of  the  agent  simulant 

Overall  damage  modes  for  chemical  submunitions  were  first  identified.  In  some  cases,  we  were 
able  to  observe  macroscopic  fracture  surfaces  that  can  indicate  fracture  initiation  sites,  crack  propagation 
direction,  types  of  loads  driving  the  cracks,  and  the  distribution  of  fragment  sizes.  At  the  end  of  this 
section,  we  discuss  the  implications  of  our  observations. 

Although  both  Type  1  and  Type  2  chemical  submunitions  were  examined,  we  concentrated  our 
efforts  on  the  damage  modes  and  fragments  in  Type  1  chemical  submunitions  (see  Figure  1-1),  which 
were  impacted  at  about  1.9  km/s.  We  examined  Type  1  submunitions  from  Test  41I-E2,  in  which  the 
target  was  impacted  at  the  center  of  mass  at  a  10-degree  angle  of  attack,  and  from  Test  41I-F4,  in  which 
the  axis  of  the  target  was  parallel  to  the  axis  of  the  projectile,  but  was  misaligned  so  that  the  projectile 
impacted  the  aft  portion  of  the  target. 

In  the  description  that  follows,  we  generally  progress  from  the  area  of  direct  impact  (i.e.,  from  the 
most  severe  damage)  to  the  area  farthest  from  the  impact  (i.e.,  to  the  least  severe  damage).  Table  2-1 
shows  some  of  the  common  overall  damage  modes.  Photographs  of  the  damaged  submunititms  appear  in 
Appendix  B  of  Mullins  et  al.  (1992).  Each  damage  mode  is  discussed  below. 

In  all  three  tests,  the  submunitions  in  the  vicinity  of  the  impact  were  heavily  fragmented  to  the 
extent  that  the  fragments  could  not  be  identified  with  the  individual  submunitions  from  which  they  came. 
Most  of  the  fragments  from  both  submunitions  were  "potato  chip-shaped"  portions  of  the  case  wall,  but 
some  fragments  were  portions  of  the  base  plate  or  portions  of  the  fuze  attachment  section  at  the 
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Table  2-1.  OVERALL  DAMAGE  MODES  OBSERVED  JN  SLED  TESTS. 


Damage  Mode 


S  ooo  nn 
OOO 
0O„ 

o^UU 

Severe  fragmentation 


Circumferential 
weld  fractures 


Lateral  crushing 


Fuse  detachment 


Failure  of  burst  diaphragm 


Figure  from 

Example  from  Submunition 

Reference* 

Threat  Simulator  No. 

B-2 

(Cannot  identify 

B-11 

submunition 
threat  simulators) 

B-15 

B-3 

120 

121 

B-10 

258 

B-16 

323 

326 

B-11 

266 

B-16  : 

■/J  330V'. 

B-11  a 

265 

B-4 

125 

135 

B-12 

298 

B-11 

B-16 

-  266  V  ' 

306 

330 

B-4 

135 

130 

B-4 

131 

B-10 

251 

261 

B-14 

311 

254 

B-10 

260 

253 

261 

Comments 


Multiple  axial  fractures 


Single  axial  fracture 


Spiral  fracture 


Usually  accompanied 
by  axial  fractures  or 
transverse  crushing 


Often  accompanied 
by  other  failure  modes 


*  S.  E.  MuHIns,  G.  B.  Booth,  and  C.  S.  Fowler,  "ERINT  HiWo-Kill  Sled  Tests  Against  Chemical  Submunition  Targets  (U),"  Technical 
Report  SEAS-SDC-3390,  Teledyne  Brown  Engineering  Report  on  Contract  No.  DASG6092-C-0036  (October  1992).  (SECRET) 
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submunition  nose.  Some  fragments  were  strips,  formed  from  axial  fractures  that  ran  down  most  or  all  of 
the  length  of  the  submunitions  that  had  length-to- width  ratios  greater  than  5.  However,  strips  woe  rarer 
than  potam  chips. 

Figure  2-1  plots  the  distributions  as  cumulative  number  of  fragments,  Ng(R),  versus  fragment 
size,  R,  where  Ng(R)  is  the  number  of  fragments  with  sizes  greater  than  R.  We  believe  that  the  lack  of 
fragments  with  sizes  less  than  1  cm  is  real,  because  there  was  no  evidence  of  the  production  of  fine 
fragments.  Furthermore,  fine  fragments  would  not  be  expected  in  such  ductile  material  at  the  loading  rates 
involved. 

In  all  tests,  at  a  distance  away  from  the  direct  impact  area,  we  also  found  large  fragments  of 
submunitions  formed  by  a  few  axial  fractures,  most  of  which  began  at  the  weld  zone  at  the  (fuze)  end 

of  the  submunition  and  ran  down  all  or  part  of  the  length  of  the  submunition.  In  some  submunitions,  axial 
fractures  did  not  propagate  the  entire  length  of  the  submunition  or  were  joined  with  spiral  fractures.  We 
believe  that  the  axial  cracks  are  iiutiated  by  ovalization  or  denting  of  the  end  of  the  cylinder. 

Another  common  damage  mode  is  a  circumferential  failure  at  the  shoulder  of  the  submunition, 
where  the  top  is  welded  to  the  cylinder.  The  top  can  be  separated  from  the  cylinder  by  axial  impact,  impact 
at  an  oblique  angle  to  the  axis,  or  lateral  impact.  In  axial  impact,  a  fragment  hits  the  fuze  axially  and 
deforms  it  plastically.  The  resulting  axial  loads  induce  bending  in  the  top  cover  and  expansion 
(mushrooming)  of  the  end  of  the  cylindo*  until  the  welded  joint  cracks  along  the  circumference  and  the 
cover  separates.  These  circumferential  fractures  apparently  occurred  in  all  the  submunitions  that  were 
either  severely  fragmented  or  that  suffered  at  least  one  axial  fracture  over  the  entire  length  of  the  cylinder. 
As  the  direction  of  impact  shifts  from  an  axial  orientation  to  an  oblique  angle  to  the  axis,  the  mechanism 
shifts  from  mushrooming  to  ovalizing  the  cylinder.  This  inward  deformation  induces  circumferential 
cracks  in  the  weldment  (and  separation  of  the  cover),  as  well  as  radial  cracks  in  the  cover  itself  (and  the 
fragmentation  of  the  cover).  Separation  of  the  top  welded  cover  was  often  accoti^anied  by  circumferential 
failure  at  the  weld  between  the  cylindrical  pcotion  of  the  main  body  and  the  bottom  of  the  submunition.  In 
a  direct  lateral  impact  (fragment  impacting  at  about  midheight  with  the  impact  velocity  along  a  radial 
direction),  the  top  can  be  removed  by  bending  along  the  axial  direction  over  the  top  piece. 

The  most  likely  mechanism  for  the  separation  of  the  bottom  closure  plate  is  the  ovalization  of  the 
submunition  case  by  a  lateral  impact  or  by  impact  at  an  oblique  angle  to  the  axis  similar  to  that  described 
above  for  removal  of  the  top  cover,  except  with  the  impact  point  near  the  bottom  of  the  case.  As  a  result  of 
ovalization,  the  brazed  ot  welded  joint  at  the  bottom  of  the  submunition  fractures  and  the  closure  plate  is 
forced  out  and/or  fractured. 

A  lateral  fragment  impact  on  a  submunition  (i.e.,  a  fragment  impacting  at  about  midheight,  with  the 
impact  velocity  along  a  radial  direction)  will  cause  ovalization  of  the  cylinder  (fiequently  accompanied  by 
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Figure  2-1 .  Fragment  size  distribution. 


circumferential  cracks  in  the  weld  and  separation  of  the  cover).  We  also  obswved  part-cylindrical  dents 
running  the  full  length  of  a  submunition,  clearly  showing  that  it  had  been  impacted  by  another  undamaged 
or  only  partially  damaged  submuiution. 

We  observed  evidence  for  several  mechanisms  that  can  lead  to  the  removal  of  the  fuze.  First,  the 
fiize  can  be  directly  impacted  by  a  fragment,  and  the  resulting  compression  and  bending  loads  force  it  out 
of  the  threaded  regicm  at  the  top  of  the  submunition.  The  fuze  then  bears  some  dents,  as  evidence  of  the 
impact,  and  its  threaded  bottom  part  is  bent,  shreaded,  or  both.  A  second  mechanism  involves  no  direct 
impact  on  the  fuze  but  an  oblique  impact  on  the  top  of  the  submunition.  The  resulting  deformation  of  the 
welded  cover  ovalizes  the  threaded  fuze-cover  attachment  region  and  squeezes  out  the  fiize.  In  this  case, 
the  fuze  does  not  bear  any  evidence  of  impact,  but  its  threaded  base  is  bent  and  ovalized  and  some  of  the 
threads  may  be  stripped.  A  third  scenario  involves  a  side  impact  of  the  submunition:  the  resulting 
translational  and/or  rotational  motion  of  the  submunition  and  inertial  forces  cause  enough  bending  at  the 
threaded  base  to  force  out  the  fuze.  In  this  case,  the  fuze  sustains  very  little  damage  except  to  the  threads. 

The  last  common  damage  mode  was  failure  of  the  burst  diaphragm.  The  burst  diaphragm  can  be 
sheared  and  ejected  by  the  buildup  of  internal  pressure  (due  to  impact  and  deformation  of  the  case),  while 
the  closure  plate  remains  attached  to  the  submunition.  This  failure  mode  is  often  accompanied  by  other 
failure  modes. 

Most  of  the  observed  shock  tubes  were  completely  flattened,  probably  by  the  internal  pressure 
buildup.  Side  impact  on  the  submunition  and  inertia  can  also  cause  the  shock  tubes  to  bend  along  the  axial 
direction  and  to  separate  from  the  fuze.  The  bottom  extremities  of  some  tubes  clearly  show  dent  marks, 
which  resulted  from  impact  with  the  internal  case  wall. 

Overall,  Chemical  Submunitions  1  and  2  undergo  the  same  material  and  structural  damage  modes, 
and  we  believe  that  the  mechanisms  driving  these  modes  are  also  the  same.  One  difference  between  the 
two  submunitions  is  that,  because  of  its  stubbier  design,  the  fuze  of  Submunition  2  does  not  fail  as 
frequently  as  that  of  Submunition  1.  The  larger  case  length-to-diameter  ratio  in  Submunition  2  may  also 
make  some  failure  modes  more  likely  than  others  and  may  require  different  critical  conditions.  For 
example,  the  smaller  diameter  of  the  case  may  make  the  top  and  bottom  region  stiffen,  hence,  higher 
impulses  or  forces  may  be  required  to  separate  the  top  and  the  bottom  from  the  cylinder. 

2.2  PROCEDURES  FOR  SINGLE  SUBMUNITION  EXPERIMENTS. 

2.2.1  Submunitions. 

We  performed  tests  on  Chemical  Submunitions  1  and  2.  At  full-scale.  Chemical  Submunition  1  is 
9  cm  (3.54  inches)  in  diameter  and  32.5  cm  (12.8  inches)  long  including  the  fuze.  Chanical  Submunition 
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2  is  6.4  cm  (2.50  inches)  in  diameter  and  32.8  cm  (12.9  inches)  long.  All  full-scale  submunitions  were 
obtained  from  Battelle  Memorial  Institute,  Huntsville,  Alabama,  and  were  fabricated  in  the  same  way  as 
the  submunitions  used  in  the  sled  tests.  We  obtained  half-scale  Chemical  Submunition  2  models  from 
Battelle  or  directly  frtmi  their  subcontracted  machine  shop  located  in  Hunstville.  The  half-scale  nKxlels 
were  used  primarily  so  that  we  could  conduct  submunition  to  submuiution  in:q)act  tests  using  the  six-inch 
air-driven  gun  described  below.  For  consistency,  we  used  half-scale  models  in  other  tests  also.  We 
obtained  the  quarter-scale  Chemical  Submunition  2  tiKxiels  from  Geno’al  Research  Corporation  (GRC)  in 
Santa  Barbara,  California.  Quarter-scale  tests  were  performed  only  to  assess  scaling  effects. 

The  bottom  of  each  submunition  has  an  end  that  tapers  down  to  an  opening  sealed  with  a  burst 
diaphragm.  Fcm’  Chemical  Submunition  1,  the  diaphragm  is  4.8  cm  (1.9  inch)  in  diameter  and  for 
Chemical  Submuniticm  2  it  is  2.2  cm  (0.875  inch)  in  diameter  (full-scale).  Each  submunition  is  fitted  with 
a  shock  tube  that  is  attached  to  the  fuze.  In  actual  operation,  the  shock  tube  is  filled  with  an  energetic 
material  that  is  initiated  by  the  fuze  and  pressurizes  the  chemical  agent,  which  then  bursts  the  diaphragm. 

In  our  impact  tests,  as  in  most  of  the  sled  tests,  the  submunitions  were  filled  to  95%  of  their  capacity  with 
water  (5%  ullage)  to  simulate  the  agent.  The  total  mass  for  Chemical  Submunition  1,  including  water,  was 
6.4  kg.  The  hysdrostatic  pressure  that  produces  diaphragm  bursting  was  repOTted  by  Battelle  to  be  27  to 
41  MPa  (40(X)  to  6000  psi)  in  Chemical  Submunition  1.  We  measured  burst  pressure  in  Chemical 
Submunition  2  to  be  41  MPa  (6000  psi).  We  measured  burst  pressure  in  the  half-scale  Chemical 
Submunition  2  to  be  90  MPa  (13,200  psi).  The  burst  pressure  in  the  quarter-scale  Chemical 
Submunition  2  was  measured  by  GRC  to  be  about  31  MPa  (4600  psi). 

FuU-scale  submunitions  are  constmcted  with  welds  that  attach  the  top  end  and  the  tapered  bottom 
to  the  cylindrical  portion  of  the  submunition.  The  top  end  of  Chemical  Submunition  1  is  conical  in  shape 
with  the  simulated  fuze  plug  threaded  into  a  hole  in  the  conical  end.  Thus,  Chemical  Submunition  1  has 
sloped  shoulders  on  the  fuze  end.  In  contrast.  Chemical  Submunition  2  has  a  flat  end  to  accommodate  the 
fuze  plug;  the  shoulder  is  slightly  rounded.  The  welds  attaching  the  bottom  section  of  each  submunition 
model  generally  were  strong  enough  not  to  fail  in  our  tests.  The  top  welds  on  both  full-scale  model 
submunitions  did,  however,  fail  in  several  tests.  The  top  welds  in  full-scale  Chemical  Submunition  1 
were  not  full-penetration  welds.  Rather,  they  began  on  the  outside  surface  of  the  submunition  and 
penetrated  not  more  than  half  the  thickness  of  the  wall.  The  welds  in  full-scale  Chemical  Submunition  2 
models  were  full-penetration  welds.  To  investigate  the  effect  of  these  top  welds  on  Chemical  Submunition 
2,  two  corresponding  full-scale,  unwelded.  Chemical  Submunition  2  nxxlels  were  specially  fabricated  at 
SRI.  We  machined  them  from  thick-walled  tubing.  On  most  of  the  half-scale  models,  the  tops  were 
welded  to  the  body  with  full-penetration  welds.  The  bodies  of  a  few  half-scale  models  were  machined 
from  solid  bar  stock  so  that  the  upper  end  of  these  models  had  no  welds.  Similarly,  the  quarter-scale 
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models  had  no  welds  but  were  machined  from  solid  rods.  The  bottom  tapered  end  of  the  quarter-scale 
models  was  threaded  into  the  cylindrical  porticHi. 

2.2.2  Loading  Techniques. 

The  majority  of  the  tests  performed  were  impact  tests.  We  impacted  individual  submunitions 
with  a  plate  or  with  another  submunition.  Preliminary  tests  were  performed  with  plates  accelerated  by 
explosive;  the  details  of  this  loading  techiuque  are  presented  in  Appendix  D. 

For  the  sizes,  shapes,  and  velocities  of  impactors  needed  to  damage  the  submuniticms  with  plate 
and  submunition-to-submunition  impact,  we  found  that  the  nwst  suitable  technique  was  to  accelerate  the 
impactor  with  a  15.2-cm-diameter  (6-inch)  air-driven  gas  gun  (Figure  2-2).  In  early  air  gun  tests,  we  used 
polyethylene  bottles  filled  with  foam  as  sabots  to  hold  the  steel  plate  impactors.  However,  for  impact 
velocities  above  1(X)  m/s,  the  gun  was  operated  at  pressures  above  0.34  MPa  (50  psi).  At  these  pressures, 
we  found  that  the  sabots  deformed,  creating  misalignment  of  the  plate  with  the  target  To  prevent  this 
problem,  we  fabricated  sabots  from  aluminum.  The  aluminum  sabots  withstand  the  maximum  pressure 
for  which  the  gun  was  designed  3.4  MPa.  Sabots  hold  the  steel  impactors  until  they  exit  the  muzzle  of  the 
gas  gun.  The  sabots  were  stripped  from  the  impactors  by  means  of  a  thick  steel  ring,  which  stopped  the 
sabot  leaving  the  plate  impactCM*  to  strike  the  target 

For  body-to-body  impacts,  the  half-scale  in:q)actor  was  lightly  bonded  to  the  front  of  the  sabot 
The  sabot  was  stripped  away  from  the  impacts  by  two  thick  steel  blocks  placed  so  the  impactor  would 
pass  between  them.  For  tests  where  the  target  was  impacted  axially,  the  target  was  placed  on  its  side  in 
alignment  with  the  axis  of  the  air  gun.  In  this  orientation,  the  ullage  volume  would  be  distributed  along  the 
top  side  of  the  model.  For  tests  where  the  target  was  impacted  on  the  side  either  by  a  plate  or  by  another 
submunition,  the  target  was  positioned  on  end  so  that  the  fuze  end  was  pointing  up.  Here,  the  ulllage 
volume  was  at  the  fuze  end.  In  preliminary  gas  gun  tests,  we  determined  the  practical  maximum  velocity 
for  a  1.85-kg  plate  to  be  about  330  m/s.  More  complete  gun  calibration  data  are  given  in  Appendix  E. 

For  impact  tests,  test  parameters  included  the  mass  of  the  impacting  plate  and  the  velocity  of 
impact  Orientation  of  the  model  was  also  an  important  parameter.  As  mentitmed  above,  models  were 
loaded  either  on  the  fiize  end  or  along  the  side.  A  few  tests  were  conducted  where  the  model  was  impacted 
at  an  angle  with  the  plate  impacting  the  shoulder  of  the  model.  For  these  tests,  the  plate  was  cut-off  along 
two  parallel  seconds,  so  that  the  plate  initially  contacted  only  the  shoulder  of  the  submunition.  Limited 
impact  tests  were  also  conducted  where  the  ullage  volume  and  location  was  a  parameter. 

Experiments  were  also  conducted  on  half-scale  models  of  Chemical  Submunition  2  canisters.  We 
chose  this  scale  so  that  complete  submunitions  could  be  launched  in  the  air-driven  gas  gun.  For  the 
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experiments,  three  different  loading  durations  were  produced  by  three  different  loading  techniques:  sheet 
explosive  tests  produced  short  duration  loads,  submunition-to-submunition  impacts  produced  intermediate 
duration  loads,  and  static  compression  tests  produced  long  duration  loads. 

For  short  duration  impulsive  loading,  a  uniform  sheet  explosive  load  was  applied  over  +40®  on  one 
side  of  the  submunition,  as  shown  in  Figure  2-3.  The  fuze  was  not  loaded.  A  layer  of  half-inch-thick 
foam  neoprene  was  placed  between  the  explosive  and  the  submunition  to  suppress  any  damage  ftom  the 
direct  shock  wave.  Thus,  any  damage  observed  was  caused  by  a  later  time  structural  response.  Sheet 
explosive  thickness  was  selected  to  provide  loads  that  bracket  canister  damage  from  slight  deformation  of 
the  loaded  surface  to  diaphragm  rupture  and  fuze  removal.  Impulses  imparted  to  the  target  submunition 
from  direct  explosive  loading  are  well  known  and  reproducible  (Lindberg  and  Colton,  1970). 

For  quasi-impulsive  loading,  the  side  of  each  target  canister  was  impacted  by  another  submunition 
at  velocities  ftom  100  to  230  m/s.  This  impact  loading  is  similar  to  the  plate  impactor  tests,  except  that  the 
impactor  is  another  submunition. 

For  quasi-static  loading,  a  load  was  applied  to  two  sides  of  a  submunition  by  compressing  it 
between  two  steel  platens,  as  shown  in  Figure  2-4.  The  steel  platens  extend  over  the  entire  length  of  the 
model  so  that  the  ends  of  the  model,  which  are  stiffer  than  the  central  portion,  carry  much  of  the  crush 
load. 

2.2.3  Measurement  Techniques. 

2.2.3.1  Velocity  Measurement.  Impact  velocity  was  measured  with  an  optical  system  installed  in 
the  barrel  of  the  gun.  Two  optical  stations  were  used,  one  10.2  cm  (4  inches)  from  the  end  of  the  barrel 
and  one  12.7  cm  (5  inches)  from  the  end.  At  each  station,  a  light  source  illuminated  the  sabot  as  it  passed 
the  measurement  point.  A  fiber  optic  light  sensor  was  embedded  in  the  light  source.  When  the  light 
reflected  off  the  sabot,  the  fiber  optic  sensor  would  convert  this  light  into  a  signal  that  was  recorded  on  a 
digital  scope.  To  enhance  this  measurement  technique,  the  sabot  was  marked  with  reflective  tape.  Two 
intervals  of  reflective  tape  each  2.54  cm  long,  were  separated  by  nonreflective  (black)  tape.  As  the  sabot 
passed  the  fiber  optic  sensor,  the  alternating  light  and  dark  surfaces  were  registered,  giving  a  reliable 
measurement  of  the  exit  velocity  of  the  sabot  and  impactor. 

2.2.3.2  High-Speed  Photography  and  Video.  In  some  tests,  submunition  response  and 
imparted  velocity  were  measured  with  a  Hycam  running  at  5(X)  frame/s.  In  some  tests,  high-speed  video 
was  used  to  record  the  response.  Framing  rates  for  the  video  were  1000  to  2000  frame/s. 


Figure  2-3.  Impulsive  ioading  of  submunitions  using  sheet  explosive. 
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Figure  2-4.  Quasi-static  loading  of  submunitions;  P  =  0.92  kbar. 
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PROCEDURES  FOR  SINGLE  SUBMUNITION  CALCULATIONS. 


2.3.1  DYNA3D. 

We  performed  finite  element  calculations  of  the  submunition  for  various  loading  conditions  to  help 
in  experiment  planning  and  interpretation,  to  investigate  responses  for  loading  conditions  not  tested,  and  to 
relate  the  external  loads  to  the  internal  loads.  Most  of  these  submunition  response  calculations  were 
performed  using  the  DYNA3D  finite  element  code  developed  at  Lawrence  Livermore  National  Laboratory 
(Whirley,  1993),  DYNA3D  is  an  explicit,  nonlinear,  three-dimensional  finite  element  code  for  analyzing 
the  large  deformation  dynamic  response  of  solids  and  structures.  Hie  equations  of  motion  are  integrated  in 
time  by  using  the  central  difference  method.  Spatial  discretization  is  achieved  with  eight-node  solid 
hexahedron  (brick)  elements. 

All  of  the  DYNA3D  calculations  were  performed  for  Chemical  Submunition  2,  using  the  full-scale 
submunition  dimensions  (6.35  cm  outer  diameter).  All  of  the  calculations  had  at  least  one  plane  of 
symmetry  through  the  axis  of  the  submunition,  which  was  used  to  reduce  the  calculation  size  by  a  factor 
of  2.  The  mesh  used  in  the  calculations  contained  4864  elements  for  the  submunition  and  fuze  and  4(X)4 
elements  for  the  internal  fluid  and  is  shown  in  Figure  2-5.  The  behavior  of  the  circumferential  weld  near 
the  shoulder  of  the  submunitions  was  not  sufficiently  characterized  to  be  incorporated  into  the  model  and 
therefore  was  neglected.  This  will  result  in  conservative  estimates  of  submunition  lethality  by  eliminating 
a  potential  weak  point  in  the  structure. 

The  fuze  cap  was  modeled,  but  the  detonator  and  shock  tube  were  neglected  since  the  only  fuze 
failure  mode  considered  in  the  calculations  was  separation  of  the  fuze  from  the  submunition  body.  The 
threshold  connection  between  the  fuze  cap  and  submunition  body  was  modeled  with  a  region  of  fuze 
material  that  had  a  stiffness  and  strength  approximately  one-half  that  of  the  submunition  steel  combined 
with  a  tied  sliding  interface  with  failure.  Separation  of  the  fuze  and  submunition  is  allowed  after  the 
interface  failure  criterion  is  exceeded.  The  failure  criterion  used  in  the  calculations  was 

(T/T)2  +  (a/N)2>l  (2.1) 

where  x  =  shear  stress,  a  =  normal  stress,  and  T  =  N  =  350  MPa.  The  values  of  T  and  N  were  chosen 
based  on  comparisons  of  the  initial  simulations  with  results  of  the  plate  impact  experiments. 

The  submunition  was  modeled  using  an  isotropic,  elastic-plastic,  hydrodynamic,  constitutive 
model  (model  number  10  in  DYNA3D)  with  a  yield  strength  of  670  MPa  and  the  plastic  hardening  stress- 
strain  behavior  specified  in  tabular  form.  The  water  was  modeled  with  the  fluid  constitutive  model  (model 
number  9  in  DYNA3D)  with  a  viscosity  of  1.0  cP  and  a  linear  pressure  volume  relationship  with  a  bulk 
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(a)  Mesh  cross  section  through 
the  symmetry  plane  showing 
internal  mesh  resolution 


(b)  Mesh  reflected  about  the 
symmetry  plane  showing 
external  canister  geometry 


Figure  2-5.  Mesh  of  the  canister  and  internal  fluid  used  in  the  DYNA3D  finite 
element  calculations. 
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modulus  of  2.07  GPa.  Modeling  fluid  response  at  a  free  surface  (such  as  that  produced  by  the  ullage)  can 
produce  difficulties  in  Lagrangian  finite  element  codes  (like  DYNA3D)  if  splashing  or  significant  flow 
occurs  producing  large  distortion  of  elements. 

To  this  problem  in  the  calculations,  we  smeared  the  ullage/water  interface,  creating  a 

region  of  water  with  porosity.  Specifically,  for  a  submunition  with  5%  ullage,  the  model  included  a 
region  of  20%  of  the  submunition  volume  that  contained  water  with  25%  porosity  and  a  region  of  80%  of 
the  volume  that  contained  water  with  no  porosity.  Rgure  2-6  shows  a  comparison  between  a  D YN A3D 
calculation  and  a  corresponding  experiment 

2.3.2  L2D. 

We  used  the  i  2n  code  to  model  end  impact  of  Chemical  Submunition  1  (Cooper,  1980).  All  such 
calculations  were  axisymmetric.  L2D  code  is  a  two-dimensional,  Lagrangian,  continuosuly  rezoned  code. 
The  advantage  of  T  7D  over  other  comparable  codes  is  in  its  rezoning  scheme.  It  uses  a  Free  Lagrange 
method  with  quadrilaterals  instead  of  triangles.  Triangular  elements  are  used  to  simplify  connectivity 
problems  that  occur  in  mesh  generation  and  in  the  rezoning  algorithms.  However,  the  code  minumzes  the 
number  of  triangles  by  converting  pairs  of  triangles  into  quadrilaterals.  This  produces  more  econOTiiical 
solutions  and  reduces  potential  numerical  problems  introduced  by  triangular  elements.  The  rezoning 
operations  consist  of  reconnecting  the  mesh,  deleting  nodes,  creating  nodes,  and  relocating  a  node  to  a 
better  position.  Mass,  momentum,  and  total  energy  are  conserved  in  the  rezoning. 

We  paformed  two  types  of  L2D  calculations  on  Chemical  Submunition  1.  First,  a  series  of 
calculations  was  performed  for  plate  inqiact  on  the  fuze.  Figure  2-7  shows  t3q)ical  results  of  one  such 
calculation  in  which  the  impact  velocity  was  200  tn/s.  The  characteristic  load  duration  was  80  ps.  The 
maximum  stress  occurs  at  the  top  of  the  cylindrical  portion  of  the  submunition  due  to  combined  bending 
and  hoop  tension  induced  by  flattening  of  the  conical  end. 

In  the  second  type  of  calculation,  the  fiize  was  replaced  by  a  plate  flush  with  the  fuze  opening. 
Pressure  loading  was  applied  to  the  entire  top  surface  of  this  submunition.  This  loading  represents  loads 
that  would  be  produced  by  a  debris  cloud  that  impacted  the  entire  end  of  the  submunition.  In  all  cases,  the 
maYimnm  stress  occuixed  at  the  top  of  the  cylindrical  portion  of  the  submunition. 

The  results  of  both  the  DYNA3D  and  the  L2D  code  calculations  are  presented  in  Section  3. 

2.4  SUBMUNITION  RESPONSE  TO  LOCALIZED  LOADING. 

We  performed  experiments  and  calculations  to  determine  the  response  of  localized  loading,  such  as 
plate  impact,  on  Chemical  Submunitions  1  and  2.  In  actual  scenarios,  the  plate  loading  represents  loading 
produced  by  an  array  of  fragments  on  the  target  submunition  that  is  partially  shielded  by  other 
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(a)  Side  view 


(b)  Cross-section 


Figure  2-6.  Comparison  of  DYNA3D  calculated  and  measured  deformation 
in  Test  13  (19-mm-thick  steel  plate  impact  at  204  m/s). 


Z  (cm) 


submunitions,  the  intact  of  the  end  of  one  submunition  into  the  side  of  the  in^acted  submuniton,  or  large 
fragments  from  parts  of  the  target  other  than  the  submunitions  or  from  the  projectile  itself.  We  present  the 
results  of  the  experiments  and  calculations  in  the  remainder  of  this  section.  Interpretation  of  the  results  is 
given  in  Section  3. 

2.4.1  Axial  impact. 

To  see  the  progression  of  damage  observed  in  one  type  of  localized  loading,  we  first  describe  a 
series  of  axial  impact  expaiments  on  Chemical  Submunition  1  at  progressively  higher  velocities.  In  all  the 
tests,  the  steel  impacmr  plate  was  8.9  cm  in  diameter  and  had  a  mass  of  1.85  kg. 

In  Test  G-10,  the  plate  velocity  was  175.2  m/s  (total  momentum  of  3.2  x  10^  dyne*s).  For  loading 
over  the  2.3-cm-diameter  fuze  end,  the  specific  impulse  was  78(X)  ktaps.*  We  used  high-speed 
photography  to  measure  a  post-impact  velocity  of  the  submunition  of  43.5  ra/s,  which  corresponds  to  a 
total  impulse  of  2.8  x  10^  dyne*s  imparted  to  the  submunition.  In  this  test,  shown  in  Figure  2-8,  no 
pieces  broke  away  from  the  body  of  the  submunition  and  all  of  the  water  remained  within  the  submunition. 
As  for  all  tests  in  tins  series,  conservation  of  momentum  shows  that  the  plate  must  have  rebounded  from 
the  submunition  on  inq)act 

The  intact  caused  sevo’e  deformation  of  the  fuze,  which  expanded  significantly  and  shortened 
from  6.4  cm  initial  length  to  4.5  cm  average  final  length  (1.9  cm  shortening),  as  shown  in  the  figure.  In 
this  test,  the  impact  also  caused  the  cone  to  flatten  and  the  top  end  of  the  cylindrical  body  to  bulge  or 
mushroom.  The  overall  shortening  of  the  conical  end  and  cylindrical  body  was  0.2  cm;  thus,  the  total 
shortening  of  the  submunition  was  2.1  cm.  The  final  diameter  of  the  mushroomed  end  was  9. 1  cm  (2.2% 
expansion).  An  incipient  crack  in  the  weld  between  the  conical  end  and  the  cylindrical  body  extended 
more  or  less  ctmtinuously  around  half  of  the  circumference  but  did  not  extend  through  the  entire  thickness. 
The  impact  also  produced  axial  cracks  that  emanated  from  the  circumferential  crack;  these  axial  cracks  were 
spaced  about  1.5  cm  apart  but  generally  extended  only  about  0.7  cm.  The  damage  of  interest  here  appears 
to  be  near  the  threshold  of  crack  initiaticm. 

In  Test  G-9,  the  measured  plate  velocity  was  198.4  m/s  (total  momentum  of  3.7  x  10“^  dyne«s). 

PCX'  loading  over  the  2.3-cm-diameter  fuze  end,  the  specific  impulse  was  8800  ktaps.  We  used  high-speed 
photography  to  measure  the  post-impact  velocity  of  53.4  m/s,  which  corresponds  to  a  total  momentum  of 
3.4  X  10^  dyne«s  imparted  to  the  submunition,  again  indicating  plate  rebound.  In  this  test,  as  shown  in 
Figure  2-9,  no  pieces  broke  away  from  the  body  of  the  submimition  and  all  of  the  water  remained  within 
the  submunition. 


*1  ktap  =  100  MPa-|is,  1  tap  =  1  dyncs/cm^. 
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Figure  2-8.  Chemical  Submunition  1  after  fuze-end  impact  of  a  1 .85-kg 
mass  at  a  velocity  of  175  m/s  (Test  G-10). 
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Figure  2-9.  Chemical  Submunition  1  after  fuze-end  impact  of  a  1 .85-kg  mass  at  a  velocity 
of  198  m/s  (Test  G-9). 


22 


The  impact  caused  significant  deformation  of  the  fuze.  It  expanded  significantly  and  shortened 
from  6.4  cm  initial  length  to  3.6  cm  average  final  length  (2.8  cm  shortening),  as  shown  in  Figure  2-9. 

The  impact  also  caused  the  conical  top  end  of  the  submunition  to  flatten  and  the  top  end  of  the  cylindrical 
body  to  bulge  or  mushroom.  The  overall  shortening  of  the  conical  end  and  cylindrical  body  was  0.3  cm; 
thus,  the  total  shortening  of  the  submunition  was  3.1  cm.  The  final  diameter  of  the  mushroomed  end  was 
9.2  cm  (3.4%  expansion).  A  crack  in  the  weld  between  the  conical  end  and  the  cylindrical  body  extended 
more  or  less  continuously  around  the  entire  circumference  but  did  not  extend  through  the  entire  thickness. 
The  impact  also  produced  axial  cracks  that  emanated  from  the  circumferential  crack;  these  axial  cracks  were 
spaced  about  1.5  cm  apart  and  generally  extended  about  1.4  cm.  Two  of  the  axial  cracks  (about  120 
apart)  extended  5  cm  axially,  then  branched  into  twin  cracks  that  continued  to  extend  as  much  as  12  cm. 

The  response  of  interest  here,  shown  in  Figure  2-9(c),  appears  to  be  slightly  above  the  threshold  of  crack 
initiation. 

In  Test  G-8,  the  measured  plate  velocity  was  233.6  m/s  (total  momentum  of  4.3  x  10^  dyne«s). 

For  loading  over  the  2.3-cm-diameter  fuze  end,  the  specific  impulse  was  10,4(X)  ktaps.  We  used  high¬ 
speed  photography  to  measure  the  post-impact  velocity  of  63.3  m/s,  which  corresponds  to  a  total 
momentum  of  4.1  x  lO"^  dyne*s  imparted  to  the  submunition.  This  momentum  is  based  on  the  entire 
submunition  and  fluid  mass  and  is  an  upper  bound,  because  the  submunition  fragmented  and  lost  its  water 
(and  hence  an  unknown  portion  of  its  mass)  within  8  ms  of  impact  but  before  the  velocity  was  measured. 
Therefore,  the  measured  velocity  corresponds  to  SOTne  smaller  unknown  mass  and  associated  momentum. 

The  impact  caused  the  fuze  to  deform  and  fail  as  shown  in  Figure  2-10.  In  this  test,  the  impact 
also  caused  the  body  of  the  submunition  to  fracture  into  several  pieces.  The  main  piece  that  broke  away 
from  the  body  of  the  submunition  was  recovered,  but  the  remaining  pieces  from  the  conical  end  could  not 
be  found.  The  fragmentation  was  initiated  by  three  major  axial  cracks  approximately  evenly  spaced  around 
the  submunition.  These  cracks  joined  the  circumferential  crack  between  the  conical  end  and  the  cylindrical 
body.  The  three  major  axial  cracks  apparently  extended  axially  for  5  cm  and  then  branched  into  multiple 
cracks.  The  cracking  along  the  weld  region  traveled  around  the  entire  submunition,  with  failure  through 
the  entire  thickness  along  only  one-third  of  the  circumference. 

In  Test  G-7,  the  plate  velocity  was  325.6  m/s  (total  momentum  of  6.0  x  10“^  dyne»s).  For  loading 
over  the  2.3-cm-diameter  fuze  end,  the  specific  impulse  was  14,500  ktaps.  We  used  high-speed 
photography  to  measure  a  post-impact  velocity  of  87. 1  m/s.  This  value  is  also  an  upper  bound,  because 
the  photographic  coverage  showed  that  within  5  ms  of  impact  the  submunition  fragmented  and  lost  its 
water  before  the  velocity  was  measured. 

The  damage  shown  in  Figure  2-1 1  is  the  greatest  response  of  interest.  The  impact  caused  severe 
deformation  of  the  fuze.  It  shortened  until  the  tip  of  the  fuze  sheared  off  around  the  internal  bore,  resulting 
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Figure  2-10.  Chemical  Submunition  1  after  fuze-end  impact  of  a 
1 .85-kg  mass  at  a  velocity  of  234  m/s  (Test  G-8). 
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in  complete  failure  of  the  fuze,  as  shown  in  the  figure.  In  this  test,  the  impact  also  caused  the  body  of  the 
submunition  to  firacture  into  multiple  pieces.  Three  of  the  pieces  that  broke  away  fiom  the  body  of  the 
submunition  were  recovered,  but  the  remaining  pieces  could  not  be  found.  The  fragmentation  was 
initiated  by  three  major  axial  cracks  approximately  evenly  ^aced  around  the  submunition.  These  cracks 
joined  the  circumferential  crack  along  the  weld  between  the  conical  end  and  the  cylindrical  body.  The  three 
major  axial  cracks  apparently  extended  axially  for  12  cm  and  then  turned  to  about  a  45®  angle  and  extended 
an  additional  3  cm. 

2.4.2  Calculations  of  Axial  Loading  Response. 

We  performed  DYNA3D  calculations  to  investigate  in  more  detail  the  response  of  Chemical 
SubmuniticMi  2  and  to  enhance  our  understanding  of  the  load-damage  reladcmship.  Because  the  load- 
damage  relaticMiship  does  not  depend  on  how  the  loads  ate  generated,  we  selected  loading  that  was  most 
convenient  to  analyze.  For  most  calculaticms,  the  load  was  specified  as  a  pressure  time  history  that  was  a 
square  pulse  (constant  pressure)  of  varying  durations. 

We  perfcMmed  one  set  of  calculations  to  investigate  the  effects  of  an  end  pressure  load  acting  on  the 
fuze  of  Chemical  Submunition  2.  The  spatial  distribution  of  the  load  was  a  flood  load  over  the  exposed 
end  area,  with  the  pressure  magnitude  scaled  by  the  cosine  of  the  surface  normal  vector  to  the  axis  of  the 
submunition.  Table  2-2  provides  a  summary  of  these  end  loading  calculations. 


Table  2-2.  Pressure  end  loading  calculations  on  Chemical  Submunition  2. 


Canister 


Calculation 

Number 

Peak 

impulse 

(ktaps) 

Load 

Duration 

(us) 

Peak 

Pressure 

(MPa) 

Shoulder 

Strain 

(®/o) 

Diaphragm 

^rain 

(%) 

CantOa 

600 

10 

600 

45 

12 

Cant Ob 

400 

10 

400 

20 

7.5 

CanlOc 

300 

10 

300 

10 

4.5 

CantOd 

2700 

900 

300 

1.3 

1.2 

CanlOe 

4050 

900 

450 

6.0 

3.0 

CanlOf 

5400 

900 

600 

35 

3.5 

CanlOg 

800 

100 

800 

13 

8.0 

Cant Oh 

1000 

100 

100 

23 

7.0 

Notes 
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The  failure  mechanism  predicted  by  these  calculations  is  a  circumferential  failure  through  the 
submunition  wall,  below  the  shoulder,  near  the  weld  location,  in  the  submunitions  used  m  the 
experiments.  The  average  plastic  strain  magnitude  through  the  submunition  wall  at  this  location  is  given  in 
column  5  of  Table  2-2.  The  critical  strain  for  failure  was  taken  as  20%. 

2.4.3  Other  Plate  Impact  Experiments. 

We  conducted  several  other  plate  impact  experiments  on  Chemical  Submunitions  1  and  2.  The 
tests  consisted  of  striking  submunitions  with  circular  steel  plates  in  prescribed  orientations  and  velocities. 
The  plates  all  had  a  diameter  of  8.9  cm.  Plate  thickness  ranged  from  9  to  38  mm  for  Chemical 
Submunition  1  and  from  6  to  24  mm  for  Chemical  Submunition  2.  Impact  orientations  included  end-on 
onto  the  fuze  end  of  the  submunition,  side-on  onto  the  center  of  the  submunition,  side-on  onto  the 
shoulder  of  the  submunition,  and  at  oblique  angles  onto  the  shoulder  of  the  submunition. 

Table  C-1  (in  Appendix  C)  summarizes  the  test  series  performed  on  full-scale  Chemical 
Submunition  1.  Most  of  the  damage  modes  observed  in  the  sled  tests  were  observed  in  the  single 
submunition  tests.  Lateral  crushing  of  the  submunition  body  was  not  observed,  since  this  mode  is 
probably  achieved  in  submunition-to-submunition  impact.  Figures  2-12  through  2-14  show  typical 
damage  to  Chemical  Submunition  1  for  three  of  the  impact  orientations  tested. 

Table  C-2  summarizes  the  test  parameters  for  the  series  of  tests  performed  on  full-scale  Qiemical 
Submunition  2.  Figures  2-15  through  2-17  show  typical  damage  to  Chemical  Submunition  2  for  three  of 
the  four  impact  orientations  considered.  Lethal  damage  was  not  demonstrated  for  end-on  impact  of  this 
submunition. 

Figures  2-18  and  2-19  summarize  the  impact  velocities  required  to  cause  fluid  leakage  in  Chemical 
Submunitions  1  and  2.  Both  submunitions  are  strongest  against  end-on  loading  because  of  the  fuze  and 
weakest  in  shoulder  loading  because  of  the  weld. 

2.4.4  Other  Plate  Impact  Calculations. 

We  performed  a  series  of  plate  impact  calculations  to  investigate  the  side  impart  response  for 
Chemical  Submunition  2  for  a  greater  variety  of  impactor  sizes,  shapes,  impact  positions,  and  velocities 
than  could  be  investigated  experimentally.  Table  2-3  gives  a  summary  of  these  calculations. 

These  calculations  were  performed  using  the  baseline  plate  dimensions,  which  had  the  same  impact 
length,  thickness,  and  mass  as  the  12-mm-thick  impact  plates  used  in  the  experiments.  The  remainder  of 
the  plate  impact  calculations  were  then  performed  with  no  uUage  volume,  which  slightly  overestimates 
damage  but  eliminates  any  variation  that  might  occur  with  ullage  volume  fraction  or  poation  as  an 


27 


Figure  2-12.  Damage  to  Chemical  Submunition  1  from  side-on  impact  at  midheight. 
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Impact  Test  21 

Impactor  Velocity  (V)  =  343  m/s 
Impactor  Mass  =  415  g 


Figure  2-14.  Damage  to  Chemical  Submunition  1  from  oblique  shoulder  impact. 
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Impact  Test  27 

Impactor  Velocity  (V)  =  200  m/s 
Impactor  Mass  =  594  g 


Impact  Orientation 


Figure  2-15.  Damage  to  welded  Chemical  Submunition  2  from  side-on  impact  at  midheight, 


Impact  Test  25 

Impactor  Velocity  (V)  =  200  m/s 
Impactor  Mass  =  598  g 


Fracture 


i'’’  '  .0^4 


Impact  Orientation 


Figure  2-16.  Damage  to  welded  Chemical  Submunition  2  from  side-on  impact  near  top 


Impact  Test  26 

Impactor  Velocity  (V)  =  205  m/s 
Impactor  Mass  =  599  g 


Figure  2-17.  Damage  to  Chemical  Submunition  2  from  oblique  shoulder  impact. 
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Figure  2-18.  Impact  velocities  to  produce  failure  in  Chemical  Submunition  1  for  different 
types  of  impacts  (1  .S-cm-thick,  8.9-cm-diameter  impactors) . 
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Figure  2-1 9.  Impact  velocities  to  produce  failure  in  Chemical  Submunition  2  for  different 
types  of  impacts  (1 .2-cm-thick,  8.9-cm-diameter  impactors). 


additional  parameter.  (All  the  other  calculations  described  in  this  report  had  the  baseline  5%  ullage 
volume.) 


The  calculations  F2,  F3,  and  F4  show  the  effect  of  impact  position  with  the  baseline  fragment 
impact.  The  impact  location  produces  approximately  a  15%  variation  in  calculated  burst  diaphragm  peak 
strain  and  a  10%  to  15%  variation  in  peak  pressure  at  the  diaphragm.  In  the  calculated  response  for  the  top 
side  impact,  the  fuze  failed;  this  did  not  occur  for  the  middle  or  bottom  side  impacts. 

2.4.5  Effects  of  Ullage. 

In  our  side  impact  experiments,  we  initially  found  a  correlation  between  the  total  volume  reduction 
and  the  failure  of  the  burst  diaphragm.  Specifically,  when  the  total  volume  reduction  was  greater  than  the 
ullage  (5%),  the  diaphragm  always  burst.  When  the  total  volume  reduction  was  less  than  the  ullage  (5%), 
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the  diaphragm  rarely  burst.  We  hypothesized  that  the  diphragm  burst  was  caused  simply  by  squeezing 
the  ullage  out  of  the  submunition  until  the  hydrostatic  pressure  could  build  up  sufficiently  to  burst  the 
diaphragm.  To  check  this  hypothesis,  we  conducted  an  experiment  on  a  full-scale  Chemical  Submunition 
2  that  had  15%  ullage.  This  experiment  (G-45)  was  performed  at  a  plate  impact  velocity  of  200  m/s,  with 
the  plate  striking  the  side  of  the  submunition  midway  along  its  length.  We  used  the  same  size  plate  and 
impact  velocity  as  that  used  in  an  earlier  experiment  (G-27)  in  which  the  rupture  diaphragm  was  burst  and 
the  submunition  had  a  volume  change  of  7.2%.  In  experiment  G-45,  the  rapture  diaphragm  was  also 
burst  and  the  volume  change  was  5.3%.  This  result  suggests  that  the  mere  squeezing  of  the  submunition 
to  a  volume  change  greater  than  the  ullage  volume  is  not  the  driving  mechanism  that  leads  to  diaphragm 
rapture.  We  hypothesize  that  it  is  the  dynamic  pressure  that  leads  to  diaphragm  rapture.  DYNA3D 
calculations  were  performed  to  check  this  hypothesis. 

The  DYNA3D  calculations  shown  in  Figure  2-20  indicate  that  the  transient  pressure  pulse  acting  on 
the  diaphragm  is  a  result  of  the  dynamic  stress  waves  in  the  internal  fluid  (water).  As  the  submunition  is 
impacted,  a  compressive  pressure  pulse  is  transmitted  through  the  canister  wall  into  the  fluid.  This 
pressure  pulse  spreads  through  the  interior  of  the  submunition,  reflecting  off  the  interior  surfaces, 
producing  a  complex  transient  pressure  history.  For  an  impact  occuring  at  midheight  on  the  submunition, 
it  takes  approximately  50  ps  for  the  pressure  pulse  to  travel  to  the  burst  diaphragm  at  the  bottom  of  the 
submunition  and  to  the  ullage  region  at  the  top  of  the  submunition.  Any  relief  of  the  pressure  pulse 
occuring  at  the  ullage  requires  about  120  |is  to  reach  the  burst  diaphragm.  For  many  of  the  submunition 
loading  conditions  considered,  the  internal  loading  on  the  diaphragm  was  sufficient  to  cause  failure  of  the 
burst  diaphragm  before  the  arrival  of  the  relief  wave  from  the  ullage  at  the  top  of  the  submunition.  The 
burst  diaphragm  would  therefore  have  failed  before  the  ullage  volume  was  removed  and  before  uniform 
compression  of  the  fluid  would  have  occurred. 

To  determine  the  effect  of  ullage  location,  we  performed  two  experiments  on  a  half-scale  Chemical 
Submunition  2.  (Use  of  half-scale  rtKxiels  is  discussed  in  Section  4.)  In  Test  G-48,  the  half-scale 
submunition  was  oriented  nose  up  so  that  the  ullage  was  on  the  opposite  end  of  the  submunition  from  the 
rapture  diaphragm.  The  0.6-cm-thick,  4.45-cm-diameter  plate  struck  the  side  of  the  model  midway  along 
its  length  at  a  velocity  of  200  m/s.  In  this  test,  the  fuze  was  removed  and  the  rapture  diaphragm  was 
burst.  In  Test  G-62,  a  half-scale  Chemical  Submunition  2  was  oriented  nose  down  so  that  the  ullage  was 
next  to  the  rapture  diaphragm.  The  test  was  identical  to  Test  G-48  in  every  other  respect.  In  Test  G-62, 
we  found  that  the  diaphragm  burst  just  as  it  had  in  Test  G-48.  We  conclude  that  the  location  of  the  ullage 
is  not  significant  in  producing  damage  to  the  rapture  diaphragm. 
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2.5  SUBMUNITION  RESPONSE  TO  SIDE  AREA  LOADING. 


From  the  sled  tests,  it  is  apparent  that  many  submunidons  are  damaged  by  loading  over  the  entire 
side  of  the  submunition.  We  performed  experiments  and  calculations  to  investigate  the  response  of 
submtmitions  to  side  area  loading.  We  designed  the  experiments  and  performed  the  calculations  so  that  the 
results  could  be  cast  in  the  form  of  the  classic  pressure-impulse  (PI)  curves  (Abrahamson  and  Lindberg, 
1971),  as  discussed  in  Section  3. 

2.5.1  Experiments. 

Fot  impulsive  short-duration  loading,  the  rupture  diaphragm  burst,  and  the  fuze  plug  was  removed 
at  an  impulse  of  64  ktaps  on  the  half-scale  model  (corresponds  to  128  ktaps  on  a  full-scale  model),  as 
shown  in  Figure  2-21.  At  an  impulse  of  47  ktaps  (94  ktaps  full-scale),  the  diaphragm  and  fuze  were 
intact 

In  submunition-to-submunition  impact  tests  (intermediate  duration  loading),  neither  the  fuzes  nor 
the  diaphragm  of  the  impacted  target  is  ruptured  at  an  impact  velocity  of  1 1 1  m/s.  The  imparted  velocity  is 
50  m/s.  At  an  impact  velocity  of  158  m/s,  the  fuze  is  removed  from  the  impacted  target  but  the  diaphragm 
is  intact.  The  imparted  velocity  is  77  m/s.  Figure  2-22  shows  views  of  the  two  submunitions  in  the  test  at 
an  impact  velocity  of  230  rn/s.  In  both,  the  fuze  is  removed  and  the  diaphragm  is  ruptured.  Volume 
change  in  the  two  submunititMis  averages  7%.  This  is  consistent  with  side  impact  tests  using  circular  plate 
impactors,  in  which  diaphragm  rupture  does  not  occur  at  impact  velocities  less  than  200  m/s. 

In  the  first  static  (long-duratitMi)  test,  the  submunition  was  crushed  uniformly  between  two  flat 
platens  until  the  diaphragm  was  ruptured  at  a  load  of  178  to  100  kN  (40,000  to  45,000  lb).  At  this  load, 
the  submunition  was  unloaded  and  the  contact  area  measured.  A  static  load  of  200  kN  (45,(X)0  lb) 
corresponds  to  a  pressure  over  the  loaded  area  of  about  130  MPa  (19,000  psi).  This  test  was  repeated 
with  custom  platens  that  only  load  the  submunition  along  its  length  between  the  hard  end  points.  The 
platen  had  rounded  ends.  Here,  the  load  to  cause  failure  is  lower  because  the  stiff  end  closures  do  not 
carry  any  of  the  crush  load.  At  the  peak  load  of  about  80  kN  (18,000  lb),  leakage  was  produced  around 
the  fuze  plug  (see  Figure  2-5).  This  load  corresponds  to  a  pressure  over  the  loaded  area  of  92  MPa 
(13,500  psi). 

2.5.2  Calculations. 

We  performed  a  series  of  calculations  to  investigate  the  side  loading  response  of  Chemical 
Submunition  2  in  more  detail.  In  these  calculations,  the  load  distribution  was  chosen  to  approximate  the 
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Figure  2-21 .  Damage  to  submunition  from  impulsive  loading 

(equivalent  full-scale  load:  impulse  =  128  ktaps; 
pressure  =  1.2  GPa). 
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Figure  2-22.  Damage  to  submunitions  from  body-to-body  impact 

(equivaient  full-scale  load:  impulse  =  269  ktaps; 
pressure  =  125  MPa). 
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loading  produced  in  the  experiments.  Therefore,  the  pressure  (and  impulse  intensity)  distribution  was 
uniform  over  the  central  portion  and  tapered  at  the  edges,  i.e.. 


Po 

0<0f 

Poexp(-0/0d)2 

0>0f 

(2.1) 


where  0f  is  the  cent^  flat  top  angle  set  to  ±17. 1°  and  Sd  is  the  decay  angle  set  to  ±27.5®.  A  typical 
calculation  of  this  type  is  shown  in  Figure  2-23.  This  calculation  showed  diaphragm  burst  and  partial 
removal  of  the  fuze.  Table  2-4  provides  a  summary  of  the  submunition  calculations  with  side  pressure 
loads  applied.  CMtical  diaphragm  strain  was  taken  to  be  in  the  range  of  13%  to  17%. 


Table  2-4.  Pressure  side  area  loading  calculations  on  Chemical  Submunition  2. 


Calculation 

Number 

Peak 

impulse 

(ktaps) 

Load 

Duration 

(us) 

Peak 

Pressure 

(MPa) 

Diaphragm 

Pressure 

(MPa) 

Diaphragm 

Strain 

(%) 

Notes 

Can9k 

80 

10 

800 

45 

11 

Can9i 

100 

10 

100 

65 

17 

Can9h 

120 

10 

120 

75 

22 

CanSb 

100 

30 

333 

20 

4.5 

•  Slight  side  denting 

Can3c 

130 

30 

433 

40 

10 

CanSd 

170 

30 

567 

70 

20 

•  Dent  and  fuze  rotation 

CanSe 

220 

30 

733 

100 

36 

•  Fuze  failure 

Can9f 

150 

100 

150 

40 

5 

Can9d 

200 

100 

200 

80 

15 

Can9e 

250 

100 

250 

130 

30 

Can9a 

300 

300 

100 

36 

7 

Can9c 

375 

300 

125 

80 

15 

Can9b 

450 

300 

150 

130 

3 

Can9g 

810 

900 

090 

40 

7 

Can9j 

990 

900 

110 

55 

15 

Fuze  Failure 

2.6  RESIDUAL  VELOCITIES. 

We  conducted  another  series  of  experiments  on  half-scale  models  of  Chemical  Submunition  2  to 
determine  the  maximum  residual  velocity  that  a  submunition  could  sustain  following  an  impact  and  still  be 
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Figure  2-23.  DYNA  3D  calculation  of  pressure  loading 
on  the  side  of  Chemical  Submunition  2. 


in  a  condition  to  function  (fuze  and  diaphragm  intact  and  no  leakage).  This  information  can  be  used  for 
post-engagement  analysis. 

Figure  2-24  shows  the  impact  orientations.  In  Test  G-54,  an  undamaged  half-scale  submunition 
filled  to  95%  of  capacity  was  impacted  at  1 1 1  m/s  by  an  empty  submunition  that  had  been  lightly  damaged 
in  a  previous  test  High-speed  motiai  pictures  of  the  test  were  taken  and  analyzed.  From  the  movies,  we 
found  that  the  residual  velocity  of  the  target  submunition  was  49.1  tn/s.  The  residual  velocity  of  the 
impactor  was  34.9  nVs.  For  this  test,  the  impacted  submunition  was  intact  and  no  fluid  was  lost  from  it 
After  considering  the  initial  mass  and  final  mass  of  each  submuniticHi,  we  found  that  total  momentum  of 
both  submunitions  was  conserved  to  within  1%.  In  a  second  test  (G-55),  the  target  submunition  was 
impacted  at  a  velocity  of  158  in/s.  The  impact  removed  the  fiize  plug  from  the  target  and  about  half  of  the 
fluid  inside  was  removed,  as  shown  in  Figure  2-22.  Residual  velocity  of  the  target  submunition  was 
76.6  m/s.  The  impactor  had  a  residual  velocity  of  55.7  m/s.  Again,  after  considering  the  initial  mass  and 
final  mass  of  each  submunition,  we  found  that  the  total  momentum  of  both  submunitions  was  conserved  to 
within  1%.  Thus,  for  submunititMis  that  are  aligned  as  in  these  tests,  the  maximum  residual  velocity  for  a 
surviving  submunition  is  in  the  range  of  50  to  75  m/s. 

We  also  performed  tests  in  which  the  axis  of  the  impacting  submunition  was  rotated  90  degrees  to 
the  axis  of  the  impacted  submunition.  At  an  impact  velocity  of  127  m/s,  both  submunitions  deformed 
around  the  impact  area  but  neither  leaked,  as  shown  in  Figure  2-25.  The  velocity  imparted  to  the  impacted 
submunition  was  71  m/s.  An  impact  velocity  of  153  m/s  caused  the  diaphragms  to  burst  in  both 
submunitions.  The  imparted  velocity  was  85  m/s. 

We  conclude  that  Chemical  Submunition  2  can  have  imparted  velocities  of  about  65  m/s  without 
leaking  for  a  range  of  impact  orientations. 
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Figure  2-24.  Body-to-body  impact  orientations. 
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Figure  2-25.  Damage  to  submunitions  from  cross-axis  body-to-body  impact. 
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SECTION  3 


PRESSURE-IMPULSE  CURVES  FOR  CHEMICAL  SUBMUNITIONS 

3.1  BACKGROUND. 

A  useful  technique  for  presenting  the  load-damage  relationship  fcx  any  structure  responding  to  a 
simple  pulse  load  is  to  characterize  the  applied  load  by  its  peak  pressure  P  and  total  impulse  I.  It  has  been 
demonstrated  for  a  wide  variety  of  stractures  that  there  are  an  infinite  number  of  combinations  of  P  and  I 
that  produce  the  same  final  damage  (Abrahamson  and  Lindberg,  1971).  For  a  particular  damage  level, 
these  combinations  of  P  and  I  are  represented  as  an  isodamage  curve  in  pressure-impulse  coordinates. 
Figure  3-1  illustrates  that  a  typical  isodamage  curve  would  appear  as  a  square  hyperbola  in  PI  space.  For 
simple  pulse  shapes  (e.g.,  exponential  or  square  wave),  the  straight  lines  shown  in  this  PI  diagram  have 
slopes  equal  to  the  characteristic  load  duration. 

The  ccxnbinations  of  peak  pressure  and  impulse  needed  to  produce  a  given  damage  level  are 
represented  by  the  isodamage  curve.  For  example,  for  an  exponential  pressure  load  with  a  characterisdc 
load  duration  of  10  ps,  a  peak  pressure  of  1000  MPa  produces  the  damage  level  specified  by  the 
isodamage  curve  shown  in  Figure  3-1.  This  same  damage  can  also  be  produced  by  different  combinations 
of  peak  pressures  and  characteristic  load  durations  such  as  1  ps  and  4200  MPa  or  100  ps  and  420  MPa.  If 
the  point  representing  the  applied  load  falls  above  and  to  the  right  of  the  isodamage  curve,  the  resulting 
damage  would  be  greater  than  that  specified  by  the  isodamage  curve.  Conversely,  points  below  and  to  the 
left  of  the  curve  indicate  that  the  resulting  damage  would  be  less  than  that  specified  by  the  isodamage 
curve. 

In  what  follows,  we  first  investigate  the  applicability  of  PI  curves  to  the  complex  pulse  shapes 
typically  produced  in  impacts.  Then,  we  develop  the  PI  curve  for  Chemical  Submunition  2  under  side  area 
loading.  We  use  this  curve  as  a  baseline  to  compare  other  PI  points  for  a  different  loaded  area,  a  different 
load  orientation,  and  a  different  submunition.  Finally,  we  develop  a  PI  curve  from  simplified  analyses 
based  on  mechanisms  that  produce  burst  diaphragm  failure  and  fuze  removal. 
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PRESSUF 


IMPULSE  (ktap*) 


*1  tap  =  1  dyne  •  s/cm^ 


Figure  3-1 .  PI  diagram  for  a  responding  structure  such  as  a  submunition 
(pressure  =  P  e"^,  impulse  =  PT). 


3.2  PRESSURE-IMPULSE  CHARACTERIZATION  FOR  CHEMICAL 

SUBMUNITION  2. 

3.2.1  External  Loads. 

Figure  3-2  shows  the  data  from  the  experiments  for  side  loading  of  Chemical  Submunition  2  at  all 
three  load  dtirations.  All  submunitions  tested  did  not  have  the  same  diaphragm  strength.  In  particular,  the 
diaphragm  strength  for  the  half-scale  models  was  91  MPa  compared  with  the  full-scale  diaphragm  strength 
of  34  MPa.  Therefoe,  we  first  normalized  those  critical  load  data  associated  with  diaphragm  failure  to  a 
standard  diaphragm  strength  of  34  MPa  as  follows. 

Calculations  of  the  side  impact  of  Chemical  Submunition  2  wctc  performed  for  a  range  of 
impulses.  Figure  3-3  shows  the  pressure  histories  at  the  diaphragm  at  three  impulse  levels  for  300  ps 
duration  loads.  Figure  3-3(d)  shows  the  calculated  peak  pressure  versus  the  loading  impulse.  The 
pressure  versus  impulse  relationship  is  highly  nonlinear.  This  curve  was  used  to  normalize  the  impulse  to 
a  diaphragm  strength  of  34  MPa.  For  example,  the  half-scale  models  tested  had  a  strength  of  91  MPa 
which,  according  to  Figure  3-3(d),  would  require  an  external  impulse  of  385  ktaps  to  produce  failure.  If 
the  same  model  had  a  diaphragm  strength  of  34  MPa,  it  would  require  an  impulse  of  295  ktaps.  Thus,  if 
the  failure  mechanism  is  diaphragm  failiue,  the  critical  impulses  for  the  91  MPa  diaphragms  were  reduced 
to  300/380  =  76%  of  their  measured  value  (for  the  same  load  duration)  so  that  they  could  be  compared 
with  the  34-MPa  strength  diaphragms  found  in  full-scale  submunitions. 

The  impulse  associated  with  the  data  points  corresponding  to  diaphragm  failure  at  a  10  ps  load 
duration  has  been  reduced  to  76%  of  the  applied  impulse  to  account  for  the  greater  diaphragm  strength  in 
the  half-scale  model  tested  and  the  calculated  combinations  of  pressure  and  impulse  that  produce  the 
damage  of  interest  Figure  3-4  shows  the  data  as  well  as  the  calculations  from  Table  2-4.  The  curve 
shown  in  the  figure  is  a  hyperbola  that  best  fits  the  data  and  the  calculated  points. 

The  damage  level  represented  by  the  experiments  is  incipient  rupture  of  the  burst  diaphragm  at  the 
short  and  intermediate  load  durations  and  leakage  around  the  fuze  plug  at  the  long  duration  load.  For  all 
the  calculations,  the  failure  mode  is  burst  diaphragm  rupture  of  a  34  MPa  strength  diaphragm.  The 
impulse  and  duration  of  the  sheet  explosive  load  is  well  known  from  previous  similar  experiments.  The 
impulse  from  the  impact  experiment  is  deduced  directly  from  measurements  of  imparted  velocity;  the  load 
duration  for  the  submunition-to-submunition  impact  is  determined  from  DYNA3D  calculations.  The 
pressure  fix)m  the  static  crush  tests  is  the  maximum  total  force  divided  by  the  maximum  contact  area.  This 
example  illustrates  how  information  from  various  sources  can  be  used  as  input  to  the  specification  of  the 
PI  curve. 
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PRESSURE  (MPa) 


10,000 


^F.D 
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+  =  Diaphragm  Failure 
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0.34  MPa 
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F  =  Fuze  Removal 
D  =  Diaphragm  Failure 


1000.0 


IMPULSE  (ktap*) 


10000.0 


*1  tap  =  1  dyne  •  s/cm  2 

Figure  3-2.  Chemical  Submunition  2  PI  data  for  various  side  loading  experiments. 
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Figure  3-3.  Calculated  effect  of  impulse  on  internal  pressure  for  side  area  loading 
of  Chemical  Submunition  2  (t  =  300  ps). 
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PRESSURE  (MPa) 


IMPULSE  (ktap*) 


1  tap  =  1  dyne  •  s/cm  2 


Figure  3-4.  PI  curve  for  external  side  load  of  Chemical  Submuition  2. 


3.2.2  Internal  Loads. 

For  each  calculation  of  pressure  loading  applied  to  the  exterior  of  Chemical  Submunition  2,  the 
response  of  the  internal  fluid  was  also  calculated.  Our  goal  was  to  find  the  combinations  of  internal  PI 
pairs  that  correspond  to  the  ex^nal  PI  pairs  that  produce  failure.  Two  internal  pressure  pulses  that 
correspond  to  external  loads  that  produce  failure  are  shown  in  Figure  3-5. 

If  we  simply  assign  the  peak  value  as  the  peak  pressure  and  the  integral  of  the  pressure  time  history 
as  the  impulse  for  such  pulses,  there  is  no  pattern  to  the  internal  PI  pairs  that  produce  failure.  The  reason 
is  that  the  classical  PI  curves  were  develq)ed  for  simple  pulse  shapes  such  as  exponential,  triangular,  or 
square.  In  many  practical  ^plications,  such  as  the  impact  of  a  submunition,  the  actual  pulse  shapes  are 
complex  due  to  wave  reverberations  in  the  impactor  and  the  target,  such  as  shown  in  Figure  3-6.  The  PI 
characterization  could  be  applied  to  these  complex  pulse  shapes  if  the  complex  pulse  could  be  represented 
as  an  equivalent  simple  pulse  with  peak  pressure  P  and  total  impulse  I  related  to  the  complex  pulse. 

To  this  end,  we  analyzed  the  response  of  a  diaphragm,  such  as  die  burst  diaphragm  in  a  chemical 
submuntion,  to  complex  pulse  shapes.  The  details  of  the  analysis  are  presented  in  Appendix  A.  The 
analysis  shows  that  complex  pulses  can  be  represented  by  equivalent  square  waves  if  the  response  time  of 
the  diaphragm  is  taken  into  account.  We  use  this  result  below  to  characterize  the  complex  pressure  pulses 
developed  in  the  submuntions,  particularly  in  the  internal  fluid,  in  terms  of  equivalent  square  pulses  with 
specified  amplitude  and  duration.  The  analysis  of  PI  curves  for  complex  pulse  shapes  (Appendix  A) 
shows  how  to  interpret  such  pulses. 

For  the  10  ps  loading  pulse,  the  pressure  shown  in  Figure  3-5(a)  has  three  distinct  pulses  that 
occur  within  the  time  from  15  to  100  |is.  Because  the  rise  on  the  first  pulse  is  abrupt,  the  time  between 
pulses  is  less  than  the  response  time  of  the  diaphragm  (about  80  ps).  The  total  impulse  associated  with 
this  pulse  is  found  by  integrating  all  three  pulses,  i.e.,  from  15  to  90  ps  or  41  ktaps.  The  equivalent  peak 
pressure  is  then  41  ktaps/75  ps  =  55  MPa.  This  load  produced  incipient  failure  of  the  diaphragm. 

For  the  3(X)  ps  duration  load,  the  pressure  shown  in  Figure  3-5(b)  is  below  the  yield  stress  of  the 
diaphragm  until  160  ps  after  the  pressure  is  applied.  This  portion  of  the  pulse  does  not  contribute  to  the 
impulse.  The  rest  of  the  pulse  does,  however.  Therefore,  the  impulse  is  taken  as  the  integral  of  the 
pressure  from  160  to  350  ps  or  96  ktaps.  The  equivalent  peak  pressure  is  1/t  =  96  ktaps/190  ps  = 

51  MPa.  This  load  was  just  short  of  producing  diaphragm  failure. 

Proceeding  in  this  manner,  we  found  the  combinations  of  internal  pressures  and  impulses  that 
produce  the  burst  diaphragm  damage  of  interest  for  each  load  duration  and  plotted  them  in  Figure  3-7.  The 
curve  indicated  is  the  best  approximation  to  the  internal  pressures  that  cause  the  burst  diaphragm  to  fail. 
This  curve  can  be  used  to  predict  burst  diaphragm  failure  from  pressure  histories  measured  in  the  fluid  in 
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(a)  T=10ms 


(b)  X  =300  ms 

CAM-4727-116A 

Figure  3-5.  Pressure  histories  calculated  inside  Chemical 
Submunition  2  for  side  area  loading. 
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Figure  3-6.  Loading  pulses  for  PI  characterization. 
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submunitions  that  are  subjected  to  side  area  loading,  in  sled  tests,  or  from  calculated  pressure  histories  as 
part  of  the  array  response.  Although  not  demonstrated  here,  it  is  likely  that  this  curve  could  also  be  used 
to  estimate  failure  for  loading  other  than  side  area  loading.  Additional  calculations  and  experiments  would 
have  to  be  performed  to  make  this  determination. 

The  critical  load  curve  for  side  loading  of  Chemical  Submunition  2  was  used  as  a  reference  for 
other  critical  load  information  determined  experimentally  and  by  calculation.  Although  con:q)lete 
information  is  not  always  available  to  define  the  PI  curve  for  other  Qrpes  of  loads  or  other  targets,  it  is 
useful  to  compare  die  PI  points  that  are  available  to  this  reference  curve  as  shown  below. 

Figure  3-8  shows  the  PI  curve  for  localized  loading  for  two  different  plate  thicknesses  on  the  side 
of  Chemical  SubmunititMi  2.  Both  the  impulse  and  pressure  required  to  produce  burst  diaphragm  rupture 
are  about  a  factor  of  2  greater  than  for  side  area  loading  because  the  loading  area  is  reduced. 

Figure  3-9  shows  the  PI  curve  for  axial  loading  of  Chemical  Submunition  2.  The  curve  shown  is 
drawn  so  that  it  is  consistent  with  the  calculated  points  for  axial  loading  over  the  fuze  and  shoulder  of  the 
submunititm.  The  failure  criterion  is  a  critical  plastic  strain  of  15%,  which  is  produced  below  the 
shoulder.  The  single  data  point  is  for  loading  on  the  end  of  the  fuze.  The  data  point  is  above  the  curve 
because,  in  the  experiment,  only  the  end  of  the  fuze  was  loaded.  In  either  case,  the  critical  loads  for  axial 
loading  exceed  those  for  side  area  loading  because  the  submunition  is  stronger  in  this  direction  and 
because  the  loading  area  is  reduced. 

3.3  PRESSURE-IMPULSE  CHARACTERIZATION  FOR  CHEMICAL 

SUBMUNITION  1. 

Figure  3-10  shows  data  for  side  area  loading  of  Chemical  Submunition  1.  The  impulsive 
asymptote  is  about  a  factor  of  3  greater  than  that  for  Chemical  Submunition  2  because  of  the  greater  wall 
thickness  of  Chemical  Submunition  1. 

Figure  3-11  shows  PI  data  for  localized  plate  loading  cm  the  side  of  Chemical  Submunition  1.  By 
comparing  these  data  with  those  of  Figure  3-8,  we  see  that  the  pressure  and  impulse  is  about  a  factor  of  2 
greater  for  Chemical  Submunition  1  than  for  Chemical  Submunition  2  because  of  the  greater  wall  thickness 
of  Chemical  Submunition  1. 

Figure  3-12  shows  PI  curves  for  axial  loading  on  Chemical  Submunition  1.  Two  types  of 
calculations  were  performed  with  the  L2D  code.  One  curve  shows  a  typical  calculation  for  impact  loading 
on  the  fuze  only,  such  as  would  cmcur  by  a  large  single  fragment.  The  failure  criterion  of  15%  effective 
plastic  strain  was  produced  below  the  shoulder  of  the  submunition.  These  calculations  agree  with  the  data 
shown.  Calculations  were  also  performed  for  loading  over  the  entire  end  of  the  submunition  in  the 
absence  of  the  fuze.  This  configuration  more  closely  resembles  an  array  of  fragments  impacting  the  end  of 
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PRESSURE  (MPa) 


IMPULSE  (ktap*) 

*1  tap  =  1  dyne  •  s/cm  2 


Figure  3-1 0.  PI  data  for  side  area  loading  of  Chemical  Submunition  1 . 
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PRESSURE  (MPa) 


IMPULSE  (ktap*) 


*1  tap  =  1  dyne  •  s/cm  2 


Figure  3-1 1 .  PI  data  for  side  impact  plate  loading  experiments 
on  Chemical  Submunition  1 . 
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the  submunition.  The  damage  mode  was  similar  to  that  for  loading  on  the  fuze  only.  The  differences  in 
area  result  in  greats  critical  loads  for  the  fuze  only  loading. 

The  critical  damage  information  presented  above  is  fOT  specific  targets  under  specific  loading 
conditions.  For  lethality  assessment,  it  is  necessary  to  predict  the  response  for  a  variety  of  targets  under  a 
range  of  loads.  This  requires  simplified  ntethods  fw  determining  the  critical  loads  for  a  range  of  targets 
and  loading  conditicms.  In  the  remainder  of  this  section,  we  present  two  simplified  analyses,  one  for  burst 
diaphragm  failure  and  one  for  fuze  plug  removal. 

3.4  PRESSURE-IMPULSE  CHARACTERIZATION  FOR  A  LONG  FLUID-FILLED 

CYLINDRICAL  SHELL. 

We  |)erf(HTiied  this  analysis  to  determine  the  conditions  under  which  an  impact  to  the  side  of  a 
fluid-filled  cylindrical  shell  representing  a  submunition  will  raise  the  fluid  pressure  to  a  level  sufficient  to 
fail  the  burst  diaphragm.  The  main  theoretical  approach  to  solving  this  problem  is  to  analyze  the  impact 
event  with  a  finite  element  code  using  a  reasonably  idealized  fragment  (e.g.,  plate-  or  rectangular-shaped 
plate)  as  discussed  in  Section  2.4.  A  useful  supplement  to  the  finite-element  analysis  is  a  classical  analysis 
in  which  the  contact  pressiue  history  is  further  idealized  to  be  a  rectangular  pulse.  The  classical  analysis  is 
especially  useful  for  lethality  assessment  if  a  formula  can  be  derived  for  the  pressure  and  impulse  of  the 
rectangular  pulse  required  to  generate  the  diaphragm  burst  pressure  in  the  fluid.  To  relate  the  presuure  and 
impulse  to  an  idealized  fragment,  we  first  approximate  the  duration  of  the  rectangular  pulse  as  the  time 
taken  to  flatten  the  side  of  the  cylinder  over  a  specified  sector  by  the  maximum  inward  deformation  divided 
by  half  the  impact  velocity.  Then,  we  assume  plastic  impact  to  approximate  the  fragment  momentum 
change  and  hence  the  applied  impulse.  The  average  contact  pressure  is  the  inq>ulse  divided  by  the 
duration,  if  we  assume  that  the  contact  area  is  constant  and  equal  to  the  final  value. 

In  Appendix  B,  we  present  the  details  of  a  classical  analysis  addressing  a  long  fluid-filled 
cylindrical  shell  subjected  to  a  rectangular  pressure  pulse  applied  over  a  longitudinal  strip  of  its  outer 
surface.  The  resulting  state  is  one  of  plane  strain.  The  case  of  a  loaded  length  shorter  than  the  cylinder 
length  is  incorporated  by  increasing  the  plane  strain  fluid  pressure  required  to  burst  the  diaphragm  by  the 
ratio  of  the  total  length  of  the  cylinder  to  the  loaded  length.  Applying  this  ratio  to  the  plane  strain  pressure 
is  equivalent  to  reducing  the  fluid  bulk  modulus  by  the  inverse  of  this  ratio. 

In  the  analysis,  the  cylindrical  shell  material  is  assumed  to  be  ligid-perfectly  plastic.  Deformation 
occurs  by  means  of  a  mechanism  formed  by  stationary  plastic  hinges,  as  shown  in  Figure  3-13.  The 
locations  of  the  stationary  hinges  are  based  on  experimental  observations,  finite  element  results,  and  the 
fluid  pressiue  being  the  diaphragm  burst  pressure  when  the  three  hinges  adjacent  to  the  loading  strip 
become  collinear.  At  each  plastic  hinge,  a  fully  plastic  moment  acts;  we  neglect  the  membrane  force 
interaction  on  the  moment-thrust  yield  curve.  To  simplify  the  kinematics,  we  use  the  geometry  of  the 
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chords  rather  than  the  arcs  between  plastic  hinges,  but  in  the  dynamics  we  assign  the  translational  and 
rotational  inertia  of  the  arcs  to  the  associated  chords. 

The  kineinatic  relaticms  and  the  equations  of  motion  are  reduced  to  a  single  differential  equation 
governing  the  rotation  of  the  chords  under  the  loading.  The  coefficients  in  this  equation  are  functions  of 
the  rotation,  so  a  mathematical  simplification  consisting  of  linearization  of  the  coefficients  is  introduced  to 
allow  an  explicit  solution  to  be  derived.  The  coefficient  of  the  fluid  pressure  buildup  is  also  linearized,  but 
the  fluid  pressure  is  held  constant  at  a  value  that  produces  the  same  compressional  strain  energy  or 
resistive  work  done. 

In  the  analysis,  formulas  are  derived  that  determine  the  pressure  and  in:q}ulse  required  to  drive  the 
three  plastic  hinges  under  the  load  into  a  collinear  position  in  which  rotational  motion  ceases.  Two  static 
collapse  pressures  arise  in  this  problem.  The  first  is  the  step  pressure  required  to  activate  the  hinge 
mechanism  and  the  second  is  the  smallest  step  pressure  that  drives  the  three  plastic  hinges  und^  the  load 
into  a  collinear  position  where  rotational  motion  ceases. 

Let  us  consider  a  long  steel  cylinder  with  the  following  specifications  characteristic  of  Chemical 
Submunition  2: 


Mean  radius 

a 

2.8575  cm  (1-1/8  ms) 

Mean  diameter 

d 

5.7150  cm  (2-1/4  ms) 

Wall  thickness 

h 

0.6350  cm  (1/4  m) 

Bilinear  yield  stress 

Oy 

680  MPa  (100,000  psi) 

Steel  density 

P 

7.80  g/cm^ 

Mass 

m  =  ph 

4.953  g/cm2 

the  initial  plastic  hinge  mechanism  as  having 

Chord  AB  angle 

®0 

65“ 

Chord  BC  angle 

<Po  =  n/2  -  00 

25“ 

Subtended  angle  of  AB 

Vo  =  Jt  -  200 

SOP 

Chord  AB  length 

ll  =  d  COS0O 

2.4153  cm 

Chord  BC  length 

I2  =  d  sin0o 

5.1795  cm 

ChOTd  AB  height 

Zl  =  ll  COS0O 

1.0207  cm 

Chord  BC  height 

Z2  =  I2  costpo 

4.6942  cm 

Combined  heights 

D  =  Zl  +  Z2  =  d 

5.7150  cm 

Breadth  of  chords 

B  =  ll  sin0o 

2.1890  cm 

For  these  values,  we  find  the  asymptote  to  be  Ps  =  94  MPa  and  Iq  =  124  ktaps. 
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Figure  3-14  con^ares  the  PI  curve  derived  from  this  analysis  with  that  obtained  from  the 
experiments  and  finite  element  analysis.  The  classical  analysis  presented  here  predicts  the  critical  load 
curve  reasonably  well  and  allows  the  PI  curve  for  diaphragm  failure  to  be  determined  readily  for  other 
parameters,  such  as  a  stronger  diaphragm  or  different  wall  thickness.  For  example,  Figure  3-14  also 
shows  the  critical  loading  curve  for  a  wall  thickness  of  0.378  cm  (1/8-inch),  one-half  the  original  wall 
thickness. 

3.5  PRESSURE-IMPULSE  CHARACTERIZATION  FOR  FUZE  PLUG  RESPONSE. 

In  the  sled  tests,  a  common  damage  mode  was  dislodging  or  removal  of  the  fuze  plug.  In  our 
experiments  on  single  submunititxis,  we  found  that  this  damage  was  produced  by  side  impact  of  the  body 
of  the  submunition.  Below,  we  analyze  this  mechanism  in  detail  to  determine  the  critical  loads. 

When  a  submunition  body  is  impacted  from  the  side,  the  inertia  of  the  fuze  causes  the  fuze  to  rotate 
about  a  point  <m  the  edge  of  the  base  of  the  fuze  as  indicated  in  Figure  3-15.  This  rotation  is  resisted  by 
the  threads.  Below,  we  determine  the  critical  loads  for  this  mechanism,  i.e.,  the  combinations  of  pressure 
and  impulse  applied  to  the  side  of  the  submunition  such  that  the  fuze  plug  rotates  a  specified  amount 

For  small  angles  of  rotation  (small  6),  the  equations  of  motion  for  the  fuze  and  the  body  of  the 
submunition  are 


where 


*  ..  h  f  * 

i9e=F,fi  -Fyi.T„ 


Fx  =  mx 


Fy  =  m  2  0  -  Ft 


P(t)  A  =  (M  +  m)  X 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


* 

l0  =  the  moment  of  inertia  of  the  fuze  about  its  center  of  mass 
G  =  angle  of  rotation 

Fx  =  horizontal  reaction  force  on  the  top  of  the  submunition  casing 

Fy  =  vertical  reaction  force  on  the  top  of  the  submunition  casing 

Ft  =  vCTtical  force  applied  by  the  threads 

To  =  tOTque  applied  by  the  threads  about  the  point  of  rotation 

h/2  =  vertical  distance  from  the  center  of  gravity  of  the  fiize  to  the  top  of  the  threads 

f/2  =  horizontal  distance  from  the  center  of  gravity  of  the  fuze  to  the  edge  of  the  fuze 
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PRESSURE  (MPa) 


IMPULSE  (ktap*) 

*1  tap  =  1  dyne  •  s/cm  2 

Figure  3-14.  Comparison  of  the  PI  curves  caiculated  for  the  anaiytical  model 
of  burst  diaphragm  failure  and  the  baseline  PI  curve 
for  Chemicai  Submunition  2. 
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=  mass  of  the  fuze 

=  mass  of  the  submunition  and  contents  except  for  the  fiize 
=  pressure  history  applied  to  the  side  of  die  submunition 
=  area  of  the  side  of  the  submunition 


Substituting  (3.2),  (3.3),  and  (3.4)  into  (3.1)  gives 


where 


Equation  (3.6)  can  be  written 


*  f 

To  =  T^-Ft2 


e  =  F  -  To’ 


where 


p.  oahA 

^  ~  2(M  +  m)  le 


To'  = 


•  -In 


*  mf^ 


le  =  le  + 

A  square  loading  pulse  has  two  regimes,  i.e. 
F  =Po’ 


F  =0 


For  initial  conditions 


0(0)  =  0 


0<t^t 


0(0)  =  0 
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Equation  (3.8)  can  be  integrated  the  first  regime  to 


0  =  (P*  -  To’)t 

e  =  (P-To')  J 

At  the  transition  between  the  two  regimes  t  =  x, 

e(t)  =  (F  -  To’)  t 
0(x)  =  (F  -  To')  x2/2 
In  the  second  regime,  for  t  >  X, 

e  =  -  To’ 

0  =  -  To’ t  +  B 

At  t  =  X,  (3.9)  and  (3.12)  give  B  =  Fx. 

In  the  second  regime, 


t<x 

t^x 


and  at  t  =  X 


-T«’t2 

0  =  -^  +  Fxt  +  C 


(F-To’)$  =  -To’$  +  Fx2  +  C 


or 


C  = 


-Fx2 


Therefore, 


(3.9) 

(3.10) 

(3.11) 

(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


70 


•  P'  • 

Fuze  rotation  stops  when  0  =  0,  which,  from  (3.12),  occurs  at  t  =  The  maximum  rotation  is 


fl 

vmax  ^  ^  T*  » 


^max  —  2 


For  a  square  pulse  shape,  the  impulse  is  given  by  I'  =  P'x  and  (3.17)  can  be  written 

rVi  U 


«max-2|^To’  P'J 


To  write  (3.18)  in  the  form  of  the  asymptotic  values  of  I'  and  F, 


x->0 


(3.17) 


(3.18) 


so  that 


fl  - 
«max  —  2  To' 

Combining  (3.18)  and  (3.19)  gives  the  nondimensional  PI  curve 


1-2 


or,  noting  that 


mhA  ,  , 
^  ”  2(M  +  m)^^ 


we  can  write  (3.20),  the  PI  curve  in  dimensional  form  as 


where  from  (3.19) 


2(M  +  m)  To] 
l2  ■  mhAP  J 


Io  =  ^^^-^V2ToemaxIe 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 
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The  PI  curves  for  Chemical  Submunitions  1  and  2  were  calculated  using  the  following  values: 


Chemical  Submunition  1 

Chemical  Submunition  2 

M 

g 

6535 

3916 

m 

g 

607 

470 

h 

cm 

2.24 

1.38 

A 

cn^ 

83.1 

65.6 

To 

dyne  cm 

3.47  X  lOlO 

2.71  X  IQlO 

®max 

radians 

0.03 

0.03 

le 

gcn^ 

3035 

1267 

Figure  3-16  compares  the  PI  curves  fcs"  fuze  removal  fw  both  submunitions  with  the  baseline 
curve  for  Chemical  Submuniticm  2.  The  value  of  To  and  Omax  fhU  in  a  range  ctMisistent  with  the  observed 
behavior  of  fuze  removal  The  particular  values  chosen  are  those  that  give  best  agreement  with  the  data. 

An  example  of  the  utility  of  this  analysis  is  that  it  can  be  used  to  determine  quickly  the  effect  of 
changing  the  target  parameters.  Figure  3-17  shows  the  effect  of  reducing  the  length  of  the  fuze  on 
Chemical  Submunition  2  while  keeping  all  other  parameters  constant  The  simplicity  of  the  fuze  analysis 
and  the  ease  of  investigating  a  number  of  parameters  makes  it  a  suitable  tool  for  lethality  assessments. 
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PRESSURE  (MPa) 


IMPULSE  (ktap*) 


*1  tap  =  1  dyne  •  s/cm2 

Figure  3-16.  PI  curves  for  fuze  removal  compared  to  baseline  curve 
for  Chemical  Submunition  2. 
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PRESSURE  (MPa) 


SECTION  4 

SCALING  EVALUATION 


The  objective  of  this  investigation  was  to  establish  whether  the  impact  response  of  complete  targets 
with  arrays  of  chemical  submuniticMis  can  be  reliably  obtained  from  subscale  tests  (Colton  et  al.,  1994). 

We  first  examined  the  results  from  tests  of  complete  targets  with  submunition  arrays  conducted  at  the 
Holloman  Sled  Test  Track  (full-scale)  and  at  the  Arnold  Engineering  Development  Center  (quarter-scale). 
The  damage  modes  in  the  full-  and  quarter-scale  tests  consisted  primarily  of  fragmentation  of  the 
submunitions,  failure  of  the  burst  diaphragm,  removal  of  the  fuze  plug,  failure  of  the  weld  between  the 
case  and  the  t(^  and  bottom  of  the  submunition  (full-scale  only),  and  firacture  of  the  casing.  In  the 
specimens  we  examined,  fracture  of  the  casing  was  accompanied  by  one  of  the  other  failure  modes, 
usually  failure  of  the  burst  diaphragm.  This  dual  failure  always  occured  in  the  1/4-scale  experiments,  and 
in  the  full-scale  experiments  with  few  exceptions. 

Examination  of  the  fragments  assumed  to  be  produced  in  the  intact  area,  in  both  the  full-  and 
quarter-scale  tests,  showed  that  the  fragment  size  distributions  essentially  scale.  To  investigate  structural 
damage  modes,  we  performed  tests  in  which  plates  at  200  m/s  impacted  full-,  half-,  and  quarter-scale 
single  submunitions  on  the  side  both  at  the  midheight  and  near  the  top  of  the  submunitions.  We  found  that 
the  fuze  failure  and  diaphragm  failure  rrKxies  and  the  overall  deformation  of  the  submunitions  were 
essentially  the  same  at  all  scales.  In  the  full-scale  submunitions  with  welds,  identical  to  those 
submunitions  used  in  full-scale  sled  tests,  weld  failure  occurred  as  it  often  did  in  the  sled  tests.  We  also 
made  and  tested  full-scale  submunitions  without  welds  and  found  that  the  responses  were  nearly  identical 
to  those  in  the  1/4-scale  (unwelded)  submunition:  no  failures  in  the  case,  failure  of  the  burst  diaphragm, 
and  partial  removal  of  the  fuze  plug. 

Figure  4-1  shows  a  comparison  of  the  damage  to  full-,  half-,  and  quarter-scale  models  impacted  on 
the  side  at  midlength.  Radial  deformation  and  diaphragm  rupture  in  the  three  naodels  are  faithfully 
replicated.  The  fuze  plug  was  removed  in  the  half-scale  model;  in  the  full-  and  quarter-scale  model,  the 
fuze  plug  was  dislodged  to  cause  a  large  leak.  Table  4-1  summarizes  the  deformations  and  volume 
changes  for  ccxnparable  full  and  scale  models.  The  slightly  lower  deformation  and  volume  change  for  the 
half-scale  shoulder  impact  experiment  is  consistent  with  a  slightly  off-center  impact  on  this  test  However, 
the  combination  of  deformation  and  volume  change  is  self-consistent  Nevertheless,  it  is  clear  that,  under 
equal  impact  conditions,  the  damage  sustained  by  each  model  size  is  very  similar. 


75 


Welded  Full-Scale  Unwelded  Full-Scale  Unwelded  Quarter-Scale 


(a)  At  Center 


Welded  Full-Scale  Unwelded  Full-Scale  Quarter— Scale 


(b)  At  Shoulder 


Figure  4-1 .  Comparison  of  damage  in  different  sizes  of  Chemical  Submuinitions  2 
subjected  to  side-on  impact. 
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Table  4- 1 .  Deformation  and  volume  change  measurements  for  Chemical  Submunitions  2 
impact  tests. 


Impact 

Location 

Full-Scale 

Welded 

Full-Scale 

Unwelded 

1/2-Scale 

Welded 

1/4-S 

Unwc 

icale 

Hded 

AR/Rq 

(%) 

AV/V^ 

(%) 

AR/Ro 

(%) 

AV/V„ 

(%) 

ARm^, 

(%) 

AV/V^ 

(%) 

ARm^, 

(%) 

AV/V^ 

(%) 

l^v 

62.5 

7.2 

56.9 

6.6 

53.1 

6.7 

57.1 

6.5 

n 

44.3 

6.2 

53.4 

5.4 

40.0 

4.3 

49.2 

5.0 

AR  =  radius  change;  R  ^  =  original  radius;  AV  =  volume  change;  V  ^  =  original  volume. 


Finally,  we  examined  in  detail  the  fractures  in  the  submunition  casings  from  the  full-scale  sled  tests 
and  the  quarto-scale  gun  tests.  We  found  that,  in  the  full-scale  submunitions,  cleavage  fractures  were 
large  and  more  or  less  straight.  Such  fractures  are  common  for  thicker-walled  structures  and  can  also  be 
induced  by  the  failure  of  die  weld  at  the  top  of  the  case.  In  contrast,  the  quarter-scale  submunitions  always 
fractured  by  ductile  void  growth,  which  results  in  jagged  tearing  axial  fractures  associated  with  large 
plastic  deformations,  and  drcumfoential  shear  fractures  in  the  regions  stiffened  by  the  closures.  In  our 
tests,  we  found  that,  although  the  measured  stress-strain  properties  for  the  full-  and  quarter-scale  materials 
were  essentially  the  same,  the  effective  fracture  resistance  of  the  quarter-scale  material  was  greater  than  that 
of  the  full-scale  material  because  of  the  differences  in  fracture  mode. 

We  ctMiclude  that  the  tmly  significant  difference  in  response  between  full-  and  quarter-scale 
submunitions  was  in  fracture  of  the  casing  including  weld  fracture.  Thus,  if  lethality  requires  the  casing  to 
be  fractimed,  the  loads  required  to  produce  this  response  will  be  greater  at  quarter-scale  than  at  full-scale 
and  lethality  estimates  from  the  scaled  tests  will  be  conservative.  If  lethality  relies  on  diaphragm  failure 
and  fuze  failure,  the  dominant  failure  modes  observed,  no  significant  scaling  differences  are  expected 
between  the  full-  and  quarter-scale  tests.  In  either  case,  the  small-scale  test  results  are  useful.  They  can  be 
used  to  investigate  response  modes  that  lead  to  lethality.  At  a  minimum,  small-scale  tests  can  be  used  K> 
determine  lower  bounds  on  the  loads  required  for  lethality.  In  many  cases,  they  can  be  used  to  quantify 
the  critical  loads. 


SECTION  5 
CONCLUSIONS 


Our  work  was  performed  primarily  to  determine  the  relationship  between  the  load  and  response  of 
individual  submunitions,  particularly  those  that  are  away  from  the  impact  area  where  the  impact  velocities 
are  hundreds  of  meters  per  second.  The  load-damage  relationship  is  a  link  between  (1)  measurements  that 
can  be  made  in  sled  tests  and  calculated  in  the  most  sophisticated  codes  available  (e.g.,  internal 
submuntion  pressure),  and  (2)  the  final  damage  to  individual  submunitions.  We  conclude  that 

•  Damage  modes  in  the  submunitions  recovered  from  sled  tests  included 

-  Fragmentation  caused  by  the  large  stresses  produced  in  the  vicinity  of  the  impact 

-  Failure  of  the  weld  between  the  top  and  the  cylindrical  casing 

-  Fracture  of  the  casing 

-  Rupture  of  the  burst  diaphragm  by  the  internal  pressure 

-  Removal  of  the  fuze  caused  by  the  impact  of  fragments  and  submunitions  against  the 
submimiton  away  from  the  impact  point 

•  Rupture  of  the  burst  diaphragm  and  removal  of  the  fuze  were  the  most  prevalent  failure 
modes. 

•  The  fuzes  in  Chemical  Submunitions  1  and  2  acted  as  good  energy  absorbers  that 
protected  the  body  of  the  submunition.  The  shoulder  was  weakest,  because  the  weld 
between  the  top  and  the  main  body  of  the  submuntion  frequently  failed.  For  side 
impact,  the  damage  mechanism  was  frequently  rupture  of  the  burst  diaphragm  or  fuze 
removal. 

•  The  failure  of  the  diaphragm  is  largely  a  dynamic  phenomena  of  the  fluid.  Therefore, 
the  amount  of  ullage  (at  least  up  to  15%)  and  its  location  were  not  important  in  damage 
to  the  ruptured  diaphragm. 

•  For  Chemical  Submunition  2  for  side  loading,  the  pressure-impulse  curve  for  internal 
pressure  could  be  related  to  the  pressure-impulse  curve  for  external  loads. 

•  For  a  range  of  pulse  durations  on  Chemical  Submunition  2,  the  external  critical  loads 
changed  by  an  order  of  magnitude  while  the  internal  loads  varied  by  only  a  factor  of  2 
because  the  response  of  the  fluid  was  moderated  by  the  strength  of  the  submunition. 

•  PI  curves  derived  from  simplified  analyses  for  side  impact  that  produces  diaphagm  or 
fuze  failure  can  be  used  to  predict  responses  for  a  range  of  parameters  of  interest  The 
critical  load  for  side  loading  is  greatly  increased  if  the  height  of  the  fuze  is  decreased. 
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•  SubmunticMis  could  survive  impacts  that  imparted  velocities  up  to  about  65  m/s.  This  is 
a  direct  input  into  a  post-engagement  analysis  to  determine  the  foo^rint  on  the  ground. 

•  Most  of  the  damage  mechanisms  are  the  same  in  the  full-,  half-,  and  quarter-scale 
submunitions,  notably  diaphragm  and  fiize  failure.  Differences  in  the  damage  to  the 
submunition  casing  result^  from  differences  in  construction  (welds  versus  no  welds— 
not  a  true  scaling  (Merence)  and  because  the  quarter-scale  material  fractured  at  greater 
loads  titan  the  fUl-scale  material.  Thus  if  fracture  of  the  submunition  casing  is  a 
significant  damage  mode  of  interest,  care  must  be  taken  in  interpreting  the  results. 

We  reccunmend  that 

•  Actional  expoiments  be  carried  out,  particularly  on  Submunition  1,  to  complete  the 
critical  load  curves  for  external  and  internal  loads  and  for  multiple  impacts. 

•  Simplified  models  be  developed  for  loading  over  a  portion  of  the  side  of  the 
submunititxi,  for  axial  impact,  and  for  internal  loads  for  loa^g  other  than  side  area 
loading. 

•  The  critical  load  curves  developed  be  integrated  into  lethality  assessment  codes  such  as 
KAPPII. 

•  Addititmal  experi^nts  be  conducted  to  refine  the  critical  survival  velocity  imparted  to  an 
impacted  submunition  by  glancing  blows  and  multiple  impacts.  Post-engagement 
analyses  of  the  results  should  also  be  conducted. 
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APPENDIX  A 

PRESSURE-IMPULSE  CURVES  FOR  COMPLEX  PULSES 


INTRODUCTION 

Pressure-impulse  (PI)  representations  of  isodamage  curves  are  well  established  for  basic  structural 
conqwnents  (such  as  beams,  plates,  and  cylindrical  shells)  subjected  to  basic  pulse  shapes,  such  as 
iH«».aliyivl  fnrms  of  sudden  or  blast  loadings,  which  include  the  rectangle,  triangle,  and  exponential.  H^ie, 
we  describe  a  method  for  obtaining  the  PI  relationship  for  a  pulse  that  is  more  general  than  the  usual  basic 
shape.  To  do  this,  we  have  chosen  probably  the  simplest  generalization,  which  consists  of  a  class  of 
pulses  with  general  variations  in  pressure  about  a  rectangular  pulse.  We  have  chosen  a  rigid,  perfectly 
plastic  circular  membrane  as  our  basic  structural  component  This  choice  sprang  from  a  requirement  to 
analyze  the  thin  diaphragm  that  forms  the  base  of  a  cylindrical  submunition  containing  fluid.  The  loading 
is  caused  by  the  rise  in  fluid  pressure  when  the  submunition  receives  an  impact  from  external  objects,  such 
as  an  adjacent  submunition. 

A  rigid-plastic  material  behavior  is  reasonable  because  we  are  considering  strains  that  are  many 
times  the  yield  strain,  and  the  failure  criterion  is  taken  as  the  biaxial  elongation.  The  initial  bending 
response  has  been  ignored  in  the  analysis.  In  fact,  for  the  practical  dimensions  chosen,  the  central 
deflection  is  less  than  the  thickness  when  membrane  forces  dtxninate  bending  resistance.  A  shortcoming 
of  the  analysis  is  probably  the  failure  criterion.  We  predict  that  the  maximum  strain  occurs  at  the  center, 
whereas  experimental  results  show  failure  at  the  supports.  The  different  behavior  is  attributed  to  a  stress 
concentration  at  the  support  caused  by  localized  bending  combined  with  the  membrane  tension.  One  may 
possibly  overcome  this  shortcoming  by  empirically  reducing  the  central  elongation  strain. 

We  have  restricted  the  analysis  to  pulses  that  cause  the  membrane  to  come  to  rest  after  the  entire 
pulse  has  been  applied.  For  pulses  that  are  variations  about  a  rectangular  pulse,  the  restriction  is  not 
severe.  For  pulses  tiiat  arc  variations  about  other  standard  shapes,  such  as  a  half  sine,  the  analysis 
requires  modification  to  account  for  cases  in  which  motion  can  cease  while  the  pulse  is  still  being  applied. 

In  the  analysis,  we  have  assumed  that  the  motion  is  dominated  by  the  first  or  fundamental  spatial 
mode.  The  analysis  can  be  modified  to  include  higher  nKxles.  Probably,  three  modes  would  provide  a 
fairly  accurate  description  of  the  mecharucs. 

Fluid-structure  interaction  was  examined,  but  is  not  included  here  because  the  fluid  against  the 
membrane  is  confined  within  the  cylindrical  walls.  Spatial  pressure  equalization  occurs  because  the  fluid 
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radial  wave  transit  time  is  very  short  relative  to  the  pulse  duration,  except  for  pulses  that  approach  an  ideal 
impulse.  For  these  extreme  cases,  we  have  chosen  not  to  include  the  effect  of  fluid-structure  interaction; 
the  consequence  is  that  for  extremely  short  pulses  our  in^ulse  values  are  low. 

The  main  result  of  the  analysis  is  that  the  PI  relationship  for  the  rectangular  pulse  is  not  affected 
strongly  by  the  pressure  variations  about  the  mean  even  when  the  variations  are  large. 

STATIC  MEMBRANE  RESPONSE 

Figure  A-1  shows  the  displacements  and  in-plane  forces  of  a  rigid,  perfectly  plastic  membrane 
subjected  to  a  static  pressure.  An  element  dr  of  radius  changes  its  length  to  ds  as  a  result  of  the  vertical 
and  radial  deflections  w(r)  and  v(r),  as  shown  in  Figure  A- 1(a).  The  corresponding  radial  and  tangential 
strains  are 

Er  =  =  [(1  +  v')^  +  w'2]l/2  - 1  (A.l) 

ee=7  (A.2) 

where  primes  have  been  introduced  to  indicate  differentiation  with  respect  to  r. 

We  postulate  that  the  membrane  forces  are  equal  and  have  the  value  corresponding  to  the  yield 
stress;  that  is. 


N0  —  Nf  —  No  —  <Toh  (A. 3) 

where  Cq  is  the  yield  stress  and  h  is  the  membrane  thickness.  This  stress  state  lies  on  the  von  Mises  yield 
ellipse 

N,2-NrNe  +  Ne2  =  No2  (A.4) 

The  loading  is  proportional,  so  we  can  represent  the  associated  flow  rule  in  the  form  of  the  Levy- 
Mises  equations 


_ ^_3JL 

Nr’-N0’“N7-2No 

where  Nf ,  N0',  and  Nz'  are  the  deviatoric  forces  defined  by 

Nr'  =  Nr-Nm  N0’  =  N0-Nm  and  Nz’  =  Nz-Nm 


(A.5) 


(A.6) 
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(a)  Deflected  membrane 


dr 


Figure  A-1 .  Elemental  displacements  and  forces. 


with  the  mean  force  defined  by 


Nm  =  5(Nr  +  Ne  +  Nz)  (A.7) 

The  stress  through  the  membrane  thickness  is  ignored  (Nz  =  0),  so 

Nr’  =  |(Nr-|Ne)  Ne’ =  | (Ne  -  5 Nr)  Nz’ =  -  ^(Nr  +  Ne)  =  -  Nm  (A.8) 
and  when  Ne  =  Nr  =  No,  we  obtain  the  deviatoric  values 

N0’  =  Nr’  =  |No  Nz'  =  -|No  =  -Nm  (A.9) 

In  the  flow  rule  (A.5),  e  is  the  equivalent  plastic  strain  defined  by 

e  =  ^  (Ei^  +  £,€0  +  602)1/2  (A.10) 

V3 

for  a  plastically  inctHnpressible  material  (Er  +  £0  +  Ez  =  0).  From  (A.9),  the  deviatoric  membrane  forces 
N0'  and  Nr  are  equal,  so  the  flow  rule  (A.5)  gives 

ee  =  er  (A.11) 

and,  from  (A.  10),  the  equivalent  plastic  strain  becomes 

£  =  2£0  =  2£r  (A.12) 

The  result  (A.1 1)  is  consistent  with  the  plastic  strain  vector  being  normal  to  the  yield  ellipse  at  the 
stress  state  Nr  =  N0.  We  need  this  result  to  determine  the  radial  displacement,  as  shown  below. 

In  our  treatment  of  the  equilibrium  equations,  we  assume  at  the  outset  that  the  displacements  and 
slopes  are  small  enough  to  allow  the  use  of  the  coordinate  r  (Lagrangian)  and  the  approximate  expressions 
for  the  two  principal  curvatures 

..  1  . 

Kr  =  -  w  K0  =  -  -  W 

By  resolving  forces  perpendicular  and  tangential  te  the  membrane  element,  we  obtain 

Nrw"  +  N0jw'  =  -po  (A.13) 
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and 


(rNr)'  =  N0 


(A.14) 


Equation  (A.14)  is  satisfied  by  Ne  =  Nr  =  N©  and  this  stress  state  simplifies  (A.13)  to 


w'  = 


No 


(A.  15) 


The  sdution  of  equaticm  (A.  15)  that  satisfies  the  conditions  of  zero  deflecticMi  at  the  boundary 
r  =  a  and  zero  slope  at  the  center  is 

w  =  |(a2-i2)  =  (A.16) 

The  value  of  the  constant  k  is  generally  small  enough  to  allow  the  radial  strain  expressitm  (A.1)  to  be 
approximated  by 

Er  =  v'  +  j  w’2  (A.  17) 

Equating  the  strains  £r  and  ee  as  required  by  (A.  12)  gives 

v'--  =  -iw'2  (A.18) 
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The  solution  of  equation  (A.18)  that  satisfies  the  conditions  of  zero  radial  displacement  at  the  center 
and  the  supptMt  (r  =  a)  is 

v  =  (A.19) 

The  displacements  (A.16)  and  (A.19)  give  the  strains 

ee  =  er  =  k2(^l  (A.20) 

the  maximum  being  at  the  center.  Setting  r  =  0  in  (A.20)  gives 

maxee  =  maxCr  =  ee(0)  =  er(0)  =  k2  (A.21) 


If  the  maximum  plastic  strain  in  a  failure  criterion  is  Cf,  then  (A.21)  with  k  given  by  (A.16)  shows  that  the 
failure  pressure  pf  is 


(A.22) 
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The  central  deflection  at  this  pressure  according  to  (A.  16)  is 

wf  =  ka  =  (A.23) 

As  an  example,  a  membrane  of  thickness  h  =  0.318  cm  (1/8  in.)  and  radius  a  =  2.54  cm  (1  in.) 
made  of  steel  with  a  yield  stress  of  Oq  =  612  MPa  (90,(X)0  psi,  Cy  =  0.3%)  has  a  maximum  strain  of 
ef  =  0.8%  when  the  pressure  is  pf  =  27.2  MPa  (4000  psi);  the  central  deflection  is  wf  =  0.226  cm 
(wf  =  0.72  h).  As  an  interesting  comparison,  the  static  collapse  pressure  of  this  disk,  when  it  is  regarded 
as  a  rigid-plastic  clamped  circular  plate,  is 

Ps  =  2.82  Oo  =  27  MPa  (3966  psi) 

If  we  are  concerned  with  large  failure  strains,  the  membrane  forces  begin  to  dominate  bending  aftw  small 
deflections;  however,  the  support  hinge  circle  will  introduce  initial  strain  concentrations.  The  plate  has  a 
collapse  mechanism  when  ps  is  reached,  whereas  the  membrane  does  not  have  a  collapse  mechanism; 
additional  pressure  is  needed  to  generate  additional  deflection  cm*  strain,  which  is  a  geometric  effect. 

DYNAMIC  MEMBRANE  RESPONSE 

For  the  dynamic  response,  we  add  the  inertial  term  m^o  the  right-hand  side  of  the  static 
equilibrium  equaticm  (A.13);  m  =  ph  is  the  membrane  mass  per  unit  surface  area,  and  the  dots  denote  time 
differentiation.  Again,  setting  Ne  =  Nf  =  Nq  leads  to 

w-  c2  ^w"  +  7  =  ^  c2  =  No/m  =  Oo/p  (A.24) 

In  this  analysis,  we  ignore  the  effect  of  the  radial  inertia;  a  comparison  of  the  vertical  and  radial 
displacement  expressions  (A.16)  and  (A.19)  shows  that  v  =  0(kw)  where  the  parameter  k  in  (A-16)  is 
small. 

If  we  introduce  the  dimensionless  quantities 

x  =  r/a  y  =  w/a  x-ct/a.  q(x)  =  p(t)a/No  (A.25) 


equation  (A.24)  can  be  put  in  the  convenient  dimensionless  form 


where 


()*=a/aT  (y=a/ax 

The  initial  c«iditi(Mis  arc  taken  as  zero  displacement  and  velocity 

y(x,0)  =  y(x,0)  =  0  ~  (A.27) 

and  the  boundary  conditions  as 

y(l,t)  =  0  (A.28) 

A  solution  ID  equation  (A.26)  with  the  initial  and  boundary  conditions  (A.27)  and  (A.28)  for  a 
pulse  loading  q(t)  can  be  obtained  rcadily  in  modal  form.  The  appropriate  modal  form  for  a  loading  q(x,x) 

is  the  Fourier-Bessel  expansion 


q(x,t)=  X  qnCt)  Jo(anx) 
1 


(A.29) 


where  the  sequence  of  numbers  a„  are  roots  of  Jo(a„)  =  0.  In  (A.29),  Jo  is  a  Bessel  function  of  order 
zero.  As  a  remindw,  the  properties  of  the  orthogonal  expansion  arc 


and 


1 

J  X  Jo(xam)  Jo(xan)dx 
o 


0  m5*n 

|jl2(an)  m  =  n 


qn('C)  =  -  f  X  Jo(xan)dx 

Jl  (ttn)  Q 


For  a  load  that  is  uniformly  distributed  over  the  membrane 


qn(t)  = 


2q(x) 
ttnJ  l(ctn) 


(A.30) 
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If  we  represent  the  deflection  yCx,!)  by  the  expansicm 


oo 


y(x,t)=  ^  yn(t)Jo(anx) 


(A.31) 


along  with  the  loading  expansion  (A.29)  in  equation  (A.26)  and  equate  the  coefficients  of  the  modes 
Jo(ocnx),  we  obtain  the  set  of  ordinary  differential  equations  for  the  time  dependent  coefficients: 

yn  +  oen2y  =  qn(x)  (A.32) 


with  qn(t)  given  by  (A.30). 

The  solutions  of  equations  (A.32)  that  satisfy  the  initial  conditions  yn(x,0)  =  yn(x,0)  =  0  are 

X  X 

oCnyfiCt)  =  sinOnX  J  qn(s)  costtnsds  -  cosOnX  J  qn(s)  sinOnSds  (A.33) 

o  o 

The  velocities  of  the  modal  amplitudes  are 

.  ^  '' 

yn(x)  =  cosOnX  J  qn(s)  costtnSds  +  sinOnX  J  qn(s)  sinUnSds  (A.34) 

o  o 

We  now  restrict  our  attenticm  to  the  simplest  class  of  pulses.  These  pulses  are  such  that  the 
membrane  is  still  in  motion  when  the  pulse  ends.  The  simplest  example  of  this  class  is  a  rectangular  pulse 
of  sufficiently  high  magnitude.  We  assume  as  an  approximation  that  the  time  f  at  which  the  membrane 

comes  to  rest  is  based  on  the  first  mode,  that  is,  f  is  given  by  yi(f )  =  0.  Let  the  pulse  duration  be  x . 
Then,  for  motion  ceasing  after  the  pulse  has  been  applied,  (A.33)  and  (A.34)  beccmie 

anyn(x)  =  Cn  sinttnX  -  Sn  cosanX  (A.35) 

yn(x)  =  CnCOSttnX  +  Sn  sinanX  (A.36) 


where 


X  X 

Cn=  f  qn(x)  cosctnXdX  Sn=  J  qn(x)  sinOnXdx  (A.37) 

o  o 
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The  duration  of  motion,  t ,  is  the  solution  of  yi(T)  -  0;  that  is. 


~  Cl 
tanaix  =  •  ^ 


(A.38) 


and  substitution  of  x  in  (A.35)  provides  an  approximation  to  the  displacement  amplitudes  yn(t)  when  the 
membrane  has  come  to  rest. 

The  dimensionless  forms  of  the  strain  equations  are 


c,  =  z'  +  5y^ 


e«  =  - 


V 

^  =  a 


If  we  again  postulate  that  ^  =  £o,  (Nr  —  N©),  we  obtain 

,  /zV  1  y’^ 

ee  =[xJ=-2  X 


(A.39) 


(A.40) 


Then,  considering  only  the  contribution  of  the  first  mode  to  the  strains,  the  spatial  integration  of  (A.40) 
gives 


ee  =^ai^  [1  -  Ji^(ai)l  yi^ 


(A.41) 


as  the  strain  at  the  center  of  the  membrane,  where  it  is  greatest.  Let  this  strain  be  the  failure  strain  Ef . 
Then,  to  produce  this  strain,  the  amplitude  of  the  first  mode  when  motion  ceases  must  be 


2ef^^ _ 

^^^^"ai[l-Ji2(ai)li/2 


(A.42) 


The  first  root  of  Jo(a)  =  0  is  ai  =  2.40,  which  gives  J i(ai)  =  0.52,  so  the  value  resulting  from  (A.42)  is 
yi(r)  =  1.14  efl/2.  When  Ef  =  0.1,  we  have  yi(r)  =  0.36,  so  wi(Tf)  is  about  one-third  of  the  membrane 

radius. 

Before  treating  a  general  pulse,  we  apply  our  analysis  to  a  rectangular  pulse  as  an  niustration  of  our 
approximate  analysis  and  to  establish  results  for  general  pulses  with  shapes  that  can  be  regarded  as 
variations  about  a  rectangular  pulse. 
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RECTANGULAR  PULSE 

Our  rectangular  pulse  is  defined  by 


P(t)  = 


0 


0<t<  t 


t<  t 


q(T)  = 


0 


0<t  <t 


T>X 


(A.43) 


where  the  dimensionless  pressure  q(t)  and  time  t  arc  related  to  p(t)  and  t  through  (A.25).  The  integrals  Cn 
and  Sn  in  (A.37),  according  to  (A.30)  and  (A,43),  become 


C _ ^ 


sinoqiT 


anJl(an)  On 


Sn  - 


2q  1  -  cosanT 


anJl(an)  ttn 


(A.44) 


Thus,  (A.35)  and  (A.44)  give  the  modal  displacement  amplitudes  yn  =  yn(^  at  time  x  in  the  form 


yn  =  — 7^^ - [sinanX  sinUnX  -  coscenX  (1  -  cosUnX)] 

ttn^JlCan) 


(A.45) 


where  the  duration  of  motion  is  assumed  to  be  determined  by  the  first  mode.  Hence,  by  (A.38)  and 
(A.44),  X  is  the  solution  of 


tanaix  =  -  Ci/Si  =  -  sinaix/(l  -  cosaix) 


or 


sinaix  =  cos(aix/2) 


cosaix  =  -  sin(aix/2) 


(A.46) 


(A.47) 


Expression  (A.47)  allows  us  to  sin^lify  the  formula  for  the  first  modal  amplitude  among  the  set  (A.45)  to 


4q 

yi= — ^ 


,  sin 

ai3ji(ai)  2 


aix 


(A.48) 


By  means  of  (A.42),  we  can  relate  the  pulse  amplitude  and  duration  to  the  failure  strain.  This  relationship  is 


_  .  aix  ai^Ji(ai) 

q  sin  -s—  = - T - — 

2  2[1  -  Ji2(ai)]i/2 


efl/2 


(A.49) 
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The  numcarical  value  of  the  coefficient  of  is  1.76. 

We  are  ctMisidering  rectangular  pulses  of  magnitudes  and  durations  that  do  not  allow  the  membrane 
to  c<Mne  to  rest  until  after  the  pulse  has  ended  (x  >  t).  The  lowest  pressure  to  ensure  this  condition  occurs 
when  the  pulse  duration  and  the  motion  are  equal  (t  =  x ). 

When  we  let  x  =  x  in  (A.46)  or  (A.47),  we  find  that  x  =  it/ai.  Substitution  of  aix  =  n  in  (A.49) 
gives  the  lowest  q  that  can  be  applied  to  the  membrane  to  produce  the  failure  strain  Cf.  With  this  minimum 
q  denoted  by  ^  (A.49)  gives 

J  1/2  ^  J.76  e,i/2  (A.50) 

^  2[1  -  Ji2(Oi)]l« 

For  values  of  q  <  qo,  the  membrane  comes  to  rest  while  the  pulse  is  still  being  applied.  A  modification  of 
the  analysis  is  required  to  treat  cases  in  which  the  membrane  comes  to  rest  when  the  maximum  strain  is  Cf 

while  the  pulse  is  still  being  applied. 

The  other  extreme  of  the  set  of  rectangular  pulses  is  the  ideal  impulse  defined  by 


/  ,(,)dx  (A.51) 

o 

in  terms  of  dimensionless  pressure  and  time.  For  a  rectangular  pulse,  the  Umiting  process  applied  to  Ci 
and  Si  in  (A.44)  gives 


Cl  =  -  Si  =  0  (A.52) 

aiJi(ai) 


and  when  it  is  applied  to  (A.46)  or  (A.47),  we  obtain  the  duration  of  motion  X  —  tc/2cli.  The  an^litude  of 
the  first  mode,  obtained  from  (A.48),  becomes 


~  2io 

ai2Ji(ai) 


(A.53) 


and  according  to  (A.49),  the  ideal  impulse  required  to  produce  the  failure  strain  Ef  is 

:  _  eJ/2 

^  [l-Ji2(ai)]l/2 


(A.54) 
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A  dimensionless  form  of  the  PI  relationship  is  readily  derived  from  the  results  of  (A.49),  (A.50), 
and  (A.54).  From  (A.49)  and  (A.54) 


sin  =1  i  =  q  X 
qo  2q 


(A.55) 


Equations  (A.50)  and  (A.54)  show  that  qo^  =  ai/2,  so  after  X  and  ^  are  set  as  follows, 

X  =  ^0  It  =  i/io 


Equation  (A.55)  gives  the  PI  relationship 


X  sin(pA)  =  1 


(A.56) 


In  terms  of  physical  quantities 


^  =  P^o  it  =  VIo 


I  =  pt 


ai^Ji(ai)  5ghgi/2 
2[1  -  Ji2(ai)]i/2  a  ^ 


aiJi(ai)  Oph 
[l-Ji2(ai)]i/2  c  f 


(A.57) 


For  a  rectangular  pulse  with  the  lowest  pressure  qo,  the  pressure  ratio  is  X  =  1  and  the 
corresponding  impulse  ratio,  according  to  (A.56),  is  given  by  sinji  =  1,  that  is,  p.  =  Jc/2;  the  impulse  for 
the  pulse  with  the  lowest  pressme  is  i  =  (Ji/2)io. 

As  an  example,  consider  a  membrane  of  radius  a  =  2.54  cm  and  thickness  h  =  0.318  cm,  so  that 
a/h  =  8,  and  a  yield  stress  of  Ob  =  612  MPa  (90,000  psi).  The  lowest  pressure  and  the  ideal  impulse 
values  are  given  by  (A.57)  as 


which,  for  a  failure  strain  of  Ef  =  0.1,  result  in  po  =  42.6  MPa  (6260  psi)  and  Iq  =  0.032  bar-s. 

The  lowest  pressure  pulse  has  a  duration  of  t  =  (a/c)x  =  (a/c)(jt/a)  =  1 18  ps  and  an  impulse  of 
I  =  0.050  bar-s. 

PULSES  OF  A  CERTAIN  CLASS 

With  the  objective  of  deriving  a  PI  relationship  for  general  pulses,  we  establish  the  methodology 
by  confining  our  attention  to  pulses  that  are  probably  among  the  least  complex  class  of  pulses.  This  class 
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consists  of  pulses  that  can  be  represented  reasonably  well  by  variations  about  a  rectangular  pulse,  as 
shown  schematically  in  Figure  A-2.  For  this  case,  the  pressure  is 

P(t)  =  p  +  P(t)  (A.58) 


in  which  p  is  the  average  pressure 

p=-  j  p(t)dt  (A.59) 

t  o 

and  p  =  p  -  p  is  the  departure  from  the  average  pressure. 

We  assume  (nice  again  that  the  pulse  does  not  allow  the  membrane  to  come  to  rest  until  the  entire 
pulse  has  been  applied.  We  also  assume  that  the  deflection  and  resulting  strains  are  adequately 
approximated  by  the  fundamental  mode.  With  these  assumptions  (A.35)  and  (A.36)  give 

«iyi(x)  =  Cl  sinait  -  Si  cosaix  (A.60) 

yi(x)  =  Ci  cosaix  + Si  sinait  (A.61) 


in  which 


X  X 

Ci=  J  qi(x)  cosaix  dx  Si=  j  qi(x)  sinaix  dx  (A.62) 

o  o 

In  (A.62),  qi(x)  is  the  coefficient  of  the  first  mode  in  the  expansion  (A.29).  We  represent  its  general 
temporal  form  by  a  Fourier  series.  Let 

-  2  - 

q(x)  =  q  +  q(x)  qi  (x)  =  — 77“  [q  +  qW]  (A.63) 

aiJi(ai) 

which  is  consistent  with  (A.S8),  and  substitute  in  the  integrals  of  (A.62).  Then 

Ci  =  C  +  C  Si  =  S+S  (A.64) 
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where 


2<|  sinan  -  ^  2q  1  -  cosan  ^^^5^ 

aiJi(ai)  ai  aiJi(ai)  ai 

o  ^  2 

£  = - f -  f  q(T)  cosaitdt  S  = - -  f  q(x)  sinaixdx  (A.66) 

aiJi(ai)  i  aiJi(ai)  ^ 

According  to  (A.61),  motion  in  the  fundamental  mode  ceases  at  a  time  x  given  by 

tanaix=-Ci/Si  =  -(C  +  C)/(S  +  S)  (A.67) 

The  duration  of  motion  is  approximately  the  same  as  that  associated  with  the  rectangular  pulse  of 
magnitude  q  if  the  integrals  C  and  S  are  much  smaller  than  C  and  S. 

Let  us  represent  in  (A.66)  the  variaticm  of  pressure  about  the  mean  by  the  variable 

7(t)  =  5(x)/q  (A.68) 

so  that  C  and  S  bectxne 

C  =  — ^ —  f  r(x)cosaixdx  S  =  — ^ —  f  r(x)sinaixdx  (A.69) 

aiJi(ai)  o  ctiJi(ai)  ^ 

Then  the  criteria  for  taking  the  duration  of  motion  as  that  associated  with  the  average  or  rectangular  pulse 
are 


~  sinaix  ^  .  .  1-cosaix  -rm 

f  r(x)  cosaixdx  « - —  f  r(x)  sinaixdx  « -  (A.70) 

6  ai  o 

The  inequalities  (A.70)  suggest  that  pulses  with  several  oscillations  about  the  mean  value  (r  =  0)  will 
fulfill  these  criteria.  We  also  note  that 


X 

j  r(x)dx  =  0 
o 


(A.71) 
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We,  therefore,  asstune  that  the  duration  is  given  by 


tanaiT  =  -  C/S 


With  this  simplification,  the  amplitude  of  the  first  mode,  given  by  (A.60)  is 


4q  .  ait  .  1 

yi=  sin-^  +  — 

ai3Ji(ai)  ^  ai 


C  cos 


ai''  .  o 


+  S  sin 


aix 


(A.72) 


(A.73) 


where  C  and  S  are  defined  by  (A.69). 

In  additicm  to  the  class  of  pulses  being  described  as  those  with  pressure  variations  about  a  constant 
value,  we  need  a  description  of  the  variation  with  time.  This  is  done  by  introducing  the  variable 


T]  =  x/t  =  t/t  (A.74) 

and  reqitiring  that  r  =  r(Ti)  be  the  same  for  all  pulses.  The  integrals  C  and  S  of  (A.69)  become 

~  2o  t”  ^  ~  2q  t"  ^  — 

C  =  — -  f  r(Ti) cosaitTidTi  S=— -  f  r(T|) sinaittidn  (A.75) 

aiJi(ai)  ^  aiJi(ai)  ^ 

or 

C  =  — - ^  S= — - Rs  (A.76) 

aiJi(ai)  aiJi(ai) 

where  Rc  and  Rs  represent  the  integrals  in  (A.75)  and  i  =  q  t  is  the  impulse. 

The  ideal  inq)ulse  causing  the  same  central  deflection  y  i  is  again  given  by  (A.53).  After  replacing 
yi  and  substituting  (A.76)  in  (A.73),  we  obtain 

io  =  ^  sin  ^  +  i  cos  ^  +  Rs  sin  (A.77) 

To  complete  the  derivation  of  the  pressure-impulse  relationship,  we  assume  that  the  effect  of  the 
pulse  with  the  lowest  pressure  can  be  represented  by  a  rectangular  pulse  of  magnitude  qo.  For  this  pulse, 
the  duration  of  motion  of  the  fundamental  mode  equals  that  of  the  pulse  (t  =  t).  We  have  already 
established  that  qo  =  (ai/2)io  for  producing  the  same  strain  ef,  so  the  general  pulse  duration  is  t  =  i^  = 
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(i/io)(flo/q)(jo/5o)  =  (MA)(2/ai).  Hence,  in  (A.77),  we  substitute  ait  =  2(mA)  to  obtain  the  PI 
relaticmship 


sin(D  +  cd(Rc  COSO)  +  Rs  sinca)  =  - 

K 


(A.78) 


where 

0)  =  mA  Ii  =  iyio  X  =  ^o  0<0)<n/2  (A.79) 

and 

Rc=  j  r (Ti)  cos2ci)ndTl  Rs=  J  r(il)  sinlcondn  (A.80) 

o  o 

For  a  specified  failure  strain,  qo  and  i©  are  given  by  (A.50)  and  (A.54).  Note  that  r(Tl)  is  the  sanw  for  all 
pulses  in  this  class.  A  convenient  form  of  representation  for  r(Ti)  is  a  Fourier  series  over  the  interval  0  ^  Tj 
<  1,  as  follows 

oo 

^(11)=  X  (an  0082^011  +  bn  sin2icnTl)  (A.81) 

n=l 

an  =  2j  r(Ti)  cos2jtnTidTi  bn  =  2j  r(Tl)  sin27tnTidTl  (A.82) 

o  o 

The  coefficient  ao  is  zero  because  of  (A.71).  Using  the  series  (A.81)  in  the  integrals  (A.80)  leads  to  the  PI 
relationship 

1  .  0)2  y  - -  sino)  =  -  0<0)<7t/2  (A.83) 

n=l  (n7t)2  -  o)2j  X 

The  departure  from  the  rectangular  pulse  is  given  by  the  series  in  (A.83)  and,  because  of  the  decay  of  the 
form  (njt)2,  we  see  that  the  higher  the  frequency  of  the  variations  the  less  the  effect  on  the  membrane 

response. 
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As  an  example,  let  r(T|)  =  ai  cos27ni  as  shown  in  Figure  A-3.  Then,  the  PI  relationship  is 

1  -  sin©  =  -  0  <  ©  <  Ji/2  (A.84) 

7c2  -  ©2  j  X 

A  line  of  constant  slope  in  the  pressure-impulse  (X„|x)  plane  is  ©  =  p/X  =  constant  If  ©  =  Jc/4,  X  =  V2/(l  - 
ai/15).  If  ai  =  1/2,  X  =  1.03^2  =  1.46  and  p  =  1.15.  A  rectangular  pulse  has  ai  =  0,  so  X  =  V2  =  1.41 

and  p  =  1 . 1 1.  Although  the  variation  varies  between  1/2  q  and  3/2  q,  the  pressure  and  impulse  required  to 
produce  a  failure  strain  are  only  3.4%  and  3.6%  high©:,  respectively,  than  a  rectangular  pulse  of 
magnitude  equal  to  die  average  pressure. 

Figure  A-4  shows  the  PI  curves  from  (A.84)  with  amplitudes  of  the  first  ©mpOTal  modes  of 
-ai  =  0,  and  1/2;  ai  =  0  represents  the  rectangular  pulse  and  ai  =  1/2  represents  a  pulse  with  a  cosine 
variaticm  with  an  amplitude  equal  to  half  the  average  pressure.  The  curves  are  close  together,  the  largest 
difference  being  in  the  low  pressure  region,  where  ©  =  7t/2  and  X  =  (l-ai/3)*^.  The  value  ai  =  1/2  gives 
X  =  1  and  hence  p  =  X©  =  (1.2)(;t/2),  which  means  that  the  mean  pressure  and  impulse  must  be  about  20% 
higher  than  the  rectangular  pulse  values  to  achieve  the  same  specified  failure  strain. 

Figure  A-5  shows  a  PI  curve  for  the  burst  diaphragm  in  Chemical  Submunition  2.  The  parameters 
used  were  ai  =  1/2,  ef  =  0.1,  Oq  =  620.7  MPa,  d  =  2.22  cm,  h  =  0.120  cm,  p  =  7.8  g/cm^. 
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PRESSURE 


*1  tap  =  1  dyne  •  s/cm2 


APPENDIX  B 


DERIVATION  OF  PRESSURE-IMPULSE  RELATIONSHIP  FOR  A  LONG  FLUID- 

FILLED  CYLINDRICAL  SHELL 

This  appendix  presents  details  of  the  analysis  described  in  Section  3.5. 

DERIVATION  OF  GOVERNING  EQUATION 

Figure  B-1  shows  the  plastic  hinge  mechanism  with  the  four  hinges  A,  B,  B',  and  C  fixed  by  the 
angle  0,  initially  at  Bq.  The  chord  lengths,  i\  and  h,  are  fixed  in  length  throughout  the  motion.  The 
figure  also  shows  the  forces,  moments,  and  coordinates  required  to  establish  the  kinematics  and  equations 
of  motion  of  chtM’ds  AB  and  BC. 

The  equations  of  translational  and  rotational  motion  of  the  chord  AB  are 


Pb  -  pli  sin0  -  Ry  =  miyi 

(B.l) 

pll  COS0  +  Hi  -  Rx  =  mix 

(B.2) 

Pb  [x-  1  j  +  Ryx  -  (Rx  +  Hi)  \  h  COS0  -  2M  =  Ii0 

(B.3) 

In  (B.l)  and  (B.2),  mi  is  the  mass  associated  with  chord  AB.  It  is  the  mass  of  the  arc  AB.  In  (B.3),  Ii  is 
the  moment  of  inertia  associated  with  chord  AB.  It  is  the  moment  of  inertia  of  the  arc  AB  about  its  center 
of  gravity.  Values  of  mi  and  Ii  are  introduced  later. 

The  corresponding  equations  of  motion  of  chord  BC  are 

pl2  sintp  +  Ry  =  m2y2 

(B.4) 

pl2  costp  -  H2  +  Rx  =  ni2x 

(B.5) 

RyX  +  (Rx  +  H2)  ^  I2  costp  —  2M  =  l2^ 

(B.6) 

In  these  equations,  the  mass  m2  and  moment  of  inertia  I2  are  determined  from  arc  BC.  They  will  be 
introduced  later. 
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The  kinematic  relations  are 


yi=Y  +  yzi 


y2  =  Y  +  zi+^Z2 


x  =  2B 


where 


(B.7) 

(B.8) 

(B.9) 

(B.IO) 


zi  =  li  COS0  Z2  =  I2  cos<p  B  =  li  sin0  =  I2  sin<p 

In  (B.7)  and  (B.8),  Y  is  the  displacement  of  the  central  hinge,  A. 

Elimination  of  the  reactions  Rx  and  Ry,  given  by  (B.2)  and  (B.4),  from  (B.l)  and  (B.5)  results  in 
miyi  +  m2y2  =  Pb  (B.  1 1) 

(mi  +  m2)x  =  Hi  -  H2  +  pD  (B.  12) 

where 

D  =  zi  +  Z2  (B.13) 

Equations  (B.l  1)  and  (B.12)  could  have  been  written  down  immediately  by  considoing  the  translational 
components  of  the  whole  system.  Equations  (B.3)  and  (B.6)  governing  rotaticxi  become 


mix  ^  +  m2y2x  -  Ii0  =  -  Pb  ^x  -  ^  J  5  P*i 

mix  ^ -  m2y2x  + 129  =  4 pih  (ziz2 -  B2)  +  (Hi  +  H2)  2M 


(B.14) 

(B.15) 


After  eliminating  the  horizontal  reactions  Hi  and  H2  among  (B.12),  (B.14),  and  (B.15),  we  obtain 

.  ...  Z1Z2  ..  BD  _  ,.  .. 

(mi  +  m2)  X  m2y2  +  Ii0  Z2  + 129  zi  = 


Pb^x  -  I  Jz2-|piD(B2-ziz2)-2MD 


(B.16) 
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At  this  stage,  we  have  (B.ll)  and  (B.16)  for  6  and  Y,  after  expressing  all  accelerations  in  terms  of  6. 
From  differentiation  of  (B.7),  (B.8),  (B.9),  and  (B.IO),  we  have 


=  -|02  +  ^e 


yi=Y-^02-|e 

y2  =  Y-zi0  2-B0-^q>2-y  ip 


9  =  ^ 
^  Z2 


9  = 


Z2 


( 


,2  N 


1  - 


B0  2  +  Z10 


dCD2 

d0 


(B.17) 


V  ^  j 

Substitution  of  (B.17)  in  (B.ll)  and  (B.16)  and  elimination  of  Y  yields 

[mT'+“^2  (’  +  +  (“»  +  “2)  2?  22  +  41, Z2  +  412  ^  ] 

[(»  -  %  Yi- 4.  ]  4  (m,  4  .2)  2,Z2  4  4 

0)2 

=  +  (B  -  b)  Z2  ]Pb-4pD(B2-ziZ2)-16MD 


f  z?  ^ 

1-  -i 


Z2 


(B.18) 


where  O)  =  0  is  the  angular  velocity  of  AB.  In  (B.  18),  the  coefficients  are  functions  of  0.  Initially, 

(0  =  0o)  the  chords  AB  and  BC  are  perpendicular,  which  means  that  initially  b2  =  ziz2,  and  so  the 
coefficient  of  the  fluid  pressure  is  initially  zero.  Also,  the  initial  value  of  D  is  the  diameter  of  the  cylinder. 
We  require  the  mechanism  to  cease  rotation  when  the  three  hinges  under  the  load  are  coUinear,  that  is, 
when 

0  =  ml.  At  this  final  angle,  we  have  z\  =  0,  Z2  =  (li  -  li)l^,  B  =  h,  and  D  =  Z2  in  the  coefficients. 

The  mass  and  moment  of  inertia  expressions  required  in  (B.18)  are  obtained  from  the  arcs  AB  and 
BC  in  the  form 


mi  =  ma  (Jt  -  20o)  m2  =  ma  (20o)  m  =  ph 
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Il  =  ma^  2  (ic/2)  -  Go 


cos^  Gq 
(it/2)  -  Go  . 


sin2  Go  ‘ 

l2  =  ma^  Go  -  — - 

“o 


(B.19) 


where  p  is  the  density  of  the  shell  material,  h  is  the  wall  thickness,  and  a  is  the  radius  of  the  shell. 

SOLUTION  OF  GOVERNING  EQUATION 
Equation  (B.18)  is  of  the  form 

^  -  A(a)y  =  B(a)Pi  -  C(a)pi  -  D(a)Mi  (B.20) 

da 

where 

y  =  Q)2  a  =  G-Go  maxa  =  a=^-Go  (0<a<a) 

P^=P/m  pi  =  p/m  Ml  =  M/m  (B.21) 

The  coefficients  A,  B,  C,  and  D  are  all  positive  for  a  in  the  range  0  <  a  <  d.  To  obtain  a  tractable 
solution,  we  approximate  the  coefficient  by  forms  that  are  linear  in  the  variable  a,  that  is,  we  let 

A(a)  =  Ao  +  Aia  B(a)  =  Bo  +  Bla  C(a)  =  Cia  D(a)  =  Do  +  Dla 

where 

Al  =  (Ai-Ao)/d  Bl  =  (Bi-Bo)/a  Cl  =  (Ci-Co)/d  Dl  =  (Di-Do)/d  (B.22) 

In  (B.22),  the  subscript  zero  indicates  the  initial  value  of  the  coefficient  (value  at  a  =  0),  while  the 
subscript  unity  indicates  the  final  value  of  the  coefficient  (value  at  a  =  d).  We  have  noted  that  C©  =  0.  To 
derive  a  solution  as  a  convenient  formula,  we  keep  the  fluid  pressure  at  the  constant  value,  which  gives  the 
same  resistive  work  done  on  the  deforming  mechanism.  This  value  is  half  the  fluid  pressure  that  exists  at 
the  end  of  the  rotational  motion,  that  is, 

1  1  [l2-(li-li 

P  =  ip  =  5K^ - - -  (B.23) 

where  K  is  the  fluid  bulk  modulus. 

An  approximate  solution  of  (B.20)  satisfying  the  initial  condition  of  y(o)  =  0  (zero  initial  velocity)  is 

y  =  yioc  +  y2a^  (B.24) 
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where 


yi  =  BoPl-DoMl 


(B.25) 


2y2  =  (AoBo  +  Bl)  pi  -  (AoDq  +  D1)  Ml  -  ^  Clpl 


(B.26) 


The  solution  (B.24)  applies  while  the  external  loading  pressure  P  is  acting.  Let  the  pulse  duration 
be  T  and  let  the  rotation  of  the  chord  AB  at  this  time  be  6L  Then,  noting  that 

yi/2  =  a)  =  d0/dt  =  dc(/dt 
we  obtain  the  value  of  a  by  integrating  (B.24).  The  integral  is 


a 


T  = 


J 


_ da _ 

(xl/2  (yj  +  y2a)i/2 


(B.27) 


which  provides  two  results  depending  on  whether  y2  is  positive  or  negative.  Before  supplying  these 
results,  let  us  examine  ranges  of  the  applied  pressure,  P. 

For  the  plastic  hinge  mechanism  to  be  activated,  die  right  hand  side  of  (B.20)  must  be  positive. 
Initially,  we  have  a  =  0,  so  activation  required  BqP  -  DqM  >  0.  This  result  gives  the  classical  static 
collapse  load 

Ps  =  DoM/Bo  (B.28) 


However,  because  of  the  combined  effect  on  the  motion  of  the  increasing  fluid  pressure  (represented  by 
the  coefficient  C)  and  the  geometry  changes,  this  collapse  load  will  not  cause  rotation.  Instead  of  (B.28), 
we  define  the  static  collapse  pressiue  as  the  pressure  required  to  just  achieve  the  rotation  a  (hinges 
coUinear).  This  pressure  acts  throughout  the  motion  and  the  rotation  ceases  when  a  =  a  (6  =  tc/2).  To 
tind  this  pressure,  set  y  (a)  =  0  in  the  soluticHi  (B.24)  to  give 

yi  +  y2a  =  0  (B.29) 


as  the  equation  for  Ps,  which  we  call  Ps  to  avoid  confusion  with  (B.28).  In  view  of  the  formulas 
(B.25)  and  (B.26),  (B.29)  gives 

Bo  +f  (AoBo  +  Bl) 


Ps  = 


(B.30) 
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as  an  explicit  fwTnula  for  Ps. 

In  the  integral  (B.27),  we  have  yi  >  0  and 


y2^0  for  P|Ps+ 

^  ^  d(AoBo  +  Bl) 


from  (B.25),  (B.26),  and  (B.29).  Hence  we  have  the  two  results 

T  =  -?Hln[(za)i/2  +  (1  +  za)i/2  ]  y2>0 

yi 

iT  -zs  yri 

1-1 1 - j  I  y2<0 

1  +  za 


T  = 


(-y2)^^ 


tan 


In  (B.32)  and  (B.33),  z  =  yoJyi.  The  angular  velocity  at  time  T  is 

©2  =  yia  ■¥  y2a2 


(B.31) 


(B-32) 

(B.33) 


(B.34) 


After  the  load  has  been  applied,  the  governing  equation  is  (B.20)  with  Pl  —  0.  For  convenience 
only,  let  a  =  a  -H  p.  Then  the  governing  equation  becomes 

^  -  [(Aq  +  Ala)  +  Aipi  y  =  -  [(D©  Dla)  Ml  -t-  j  Clplfi]  -  (DlM*  +  ^  Clpl)P  (B.35) 
dp 

and  the  initial  condition  is 

y  (o)  =  ©2  (B.36) 

An  approximate  solution  is 

y  =  ©2+  yiP  +  y2p2  (B.37) 

where 

yi  =  (Ao  +  Ald)  ©2-|^(Do  +  Dla)  -h  ^  C  ^p^  a  j  (B.38) 
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(B.39) 


2y2  =  (Ao  +  Aia)|(Ao  + Ala)  -  [(Do  +  Dla)  Ml  +  |c  ipi  a  j  | 

+  Al52-JblMl  +  |c  iplj 

We  want  the  angular  velocity  to  be  zero  when  0  =  n/2,  a  =  a  =  (ic/2)  —  6o  or  P  =  p  =  a  —  a.  This 
condidcHi  is  obtained  from  the  solution  (B.37)  in  the  form 

+  yiP  +  y2P^  =  0  (B.40) 

Substituticm  in  (B.40)  of  the  expressions  (B.38)  and  (B.39)  leads  to 

|l  +  (Ao  +  Ala)p  +  ^[(Ao  +  Aia)2  +  Al]  .joSP- 

=  [l  +  ^(Ao  + Ala  )P  j[(Do  +  Dia)^Ps  +  5 Clap  ]P+5[d^^Ps  +  ^dp  jp2 

(B.41) 

where  mo^  is  given  by  (B.34).  Equation  (B.34)  with  yi  and  y2  given  by  (B.25)  and  (B.26)  is 

m©2  =  [Bo  +  5(AoBo  +  Bl)  a  ]aP-[Bo  +  ^(AoDo  +  Dl)  a  ]aPs-|;Cla2p  (B.42) 

so  (B.41)  and  (B.42)  determine  P  explicitly  for  each  choice  of  the  rotation  angle  a  in  the  range 
0  <  a  <  p.  The  resulting  value  of  P  is  compared  with  (B.31)  to  see  if  y2  >  0  or  y2  <  0  in  order  to  select 
the  appropriate  pulse  duration  formula,  (B.32)  or  (B.33). 

PRESSURE-IMPULSE  RELATIONSHIP 

The  inqiulse  of  a  rectangular  pulse  is  simply  I  =  PT.  Hence,  (B.32)  and  (B.33)  give 


y2>0 

(B.43) 

2ml/2p  A  -za  ll^ 

(-my2)l^  [1  +za 

y2<0 

(B.44) 
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The  angle  a  can  be  regarded  as  a  parameter.  After  eliminating  mS^  between  (B,41)  and  (B,42)  we  obtain 
P  =  P(a)  and  hence,  according  to  (B.25)  and  (B.26),  we  have  yi(a)  and  y2(S)>  giving  z  =  z(&).  Again, 
the  inequality  (B.31)  determines  the  choice  of  (B.43)  or  (B.44)  for  the  impulse. 

If  we  cwisider  a  sequence  of  pulses  in  which  P  — > «»  and  T  — » 0  in  such  a  way  that  I  =  PT 
decreases  but  I  I©,  an  ideal  impulse,  we  can  determine  the  value  of  lo  from  (B.43).  Because  z  tends  to 
a  constant  and  a  -»  0,  the  limiting  process  gives 


2P  ,  2Pal/2 

3w(a«)W  -  -jfTj- 


yi  Bo  ' 


(B.45) 


where  m©2  is  (B.41)  with  a  =  0  (p  =  a),  that  is, 

1^1  +  Aoa>  +  Al)  d2  jm©2  ~  +  ^Aod^Do  +  j d^Ps pd2(B.46) 


SUBMUNITION  RESPONSE 

The  mass  and  naoments  of  inertia  representing  the  arcs  are 


Chord  AB 

ml  =  a  (Jt  -  20o)m 

2.4936  m  (cm) 

masubmunition 

Chord  BC 

m2  =  a  (20o)m 

6.4835  m  (cm) 

masubmunition 

Ch(xd  AB  inertia 

^  .a/"  sin29o  ^ 

Il-2a3  9o- 

1  <Po  j 

1.2585  m  (cm3) 

Chord  BC  inertia 

m 

l2  =  2a3[eo-";'®»  ] 

19.1565  m  (cm3) 

m 

The  final  values  of  the  geometrical  quantities  are 

e  =  7t/2  9  =  sin-1  =  27.8°  \j/  =  sin’l  Ol/a)  =  57.7° 

zi  =  0  Z2  =  i2  COS9  =  4.58 17  cm  D  =  Z2  =  4.58 17  cm 

B  =  li  =  I2  sin  9  =  2.4156  cm 
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The  initial  values  of  the  coefficients  in  (B.18)  are 


dco^ 

Coefficient  of 

dO 

144.58  m 

(cm^) 

Coefficient  of  co2 

-245.58  m 

(cm^) 

Coefficient  of  P 

79.38 

cm^ 

Coefficient  of  p 

0 

Coefficient  of  M 

-91.44 

cm 

where  m  is  e)q)iessed  in  grams  per  square  centimeter.  (P  and  p  in  dyncs/crtP’,  M  in  dyne-cm  per 
centimeter,  m  in  radians  per  second). 

The  final  values  of  the  coefficients  in  (B.  1 8)  are 


d(fl2 

Coefficient  of 

dO 

71.34  m 

(cm^) 

Cloefficient  of  (fl^ 

-  50.77  m 

(cm4) 

Coefficient  of  P 

79.18 

cm3 

Coefficient  of  p 

-  106.90 

cm3 

(Coefficient  of  M 

-  73.31 

cm 

The  initial  and  final  values  of  the  coefficients  in  (B.18)  give  the  following  values  for  the 
coefficients  in  (B.20): 

Ao=1.70  Al=-2.25  Bo  =  0.55cm-1  Bl  =  1.28cm-1 

Co  =  0  Cl  =  3.44cm-1  Do  =  0.63cm-3  Dl=0.92cm-3 

and  the  final  value  of  a  is  a  =  (7t/2)  -  0o  =  0.4363. 

The  plastic  hinge  moment  is 

M  =  Oyh7/4  =  0.69  X  10^  dyne-cm/cm  (69  MPa-cm^) 

so  the  static  collapse  pressure  that  triggers  the  mechanism  is 

Ps  =  (Doffio)M  =  79MPa 

The  values  of  yi  and  y2  given  by  (B.25)  and  (B.26)  are 

myi  =0.55  (P-Ps) 
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my2  =  1.11  P- 1.15  Ps- 0.86  p 
in  which  the  maximum  fluid  pressure,  according  to  (B.23),  is 

p  =  63  MPa 

For  the  choice  of  integral,  (B.32)  or  (B.33),  we  have 

y2<0  79MPa<P<lllMPa 

y2>0  P>  111  MPa 

The  lowest  pressure  that  achieves  the  rotation  6  =  n/2  (a  =  a)  given  by  yi  +  ay2  =  0,  (or  z  =  y2/yi  =  - 
1/a  =  -2.29),  is 


Ps  =  94  MPa 

This  value  is  one  of  the  asymptotes  in  the  PI  diagram.  The  other  asymptote  given  by  (B.45)  and  (B.46)  is 

Iq  =  124  ktaps  (0.124  bar»s)  (1.24  MPa*s) 

Figure  B-2  shows  the  PI  diagram  for  this  example.  It  is  the  relationship  between  the  pressure  and 
impulse  fw  all  rectangular  pulses  applied  to  a  strip  of  a  specific  cylindrical  submunition,  when  the  pulses 
cause  the  same  maximum  pressure  in  the  contained  fluid. 


LOAD  LENGTH  EFFECT 

When  the  cylinder  is  loaded  on  a  part  of  its  length,  the  fluid  pressure  that  has  to  be  generated  in  the 
foregoing  analysis  fcM-  a  plane  strain  state  must  increased  to  achieve  the  diaphragm  burst  pressure.  To 
approximate  the  larger  pressure  required,  we  have  adopted  a  quasi-static  approach.  We  multiply  the  burst 
pressure  pd  by  the  ratio  of  two  lengths,  consisting  of  the  loaded  length  plus  the  length  between  the  loading 
and  the  diaphragm  divided  by  the  loaded  length.  If  the  loaded  length  is  H  and  the  unloaded  length  above 
the  diaphragm  is  Hd,  the  new  plane  strain  pressure  required  is 


P  = 


H  +  Hd  = 

“H  W 


(B.47) 


Essentially,  we  have  a  plane  strain  fluid  with  a  lower  bulk  modulus  equal  to 

H 


Kd  = 


H  +  Hd 


K 


(B.48) 
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PRESSURE  (MPa) 


(B.49) 


To  obtain  the  higher  pressure  p  of  (B.47),  we  need  a  larger  fluid  volumetric  strain,  £d,  given  by 


ed  = 


H  +  Hd 
H 


e 


If  we  keep  the  condition  that  the  three  plastic  hinges  are  collinear  when  the  maximum  fluid  pressure 
is  reached,  the  requirement  of  an  increased  volumetric  strain  introduces  a  small  ch^ge  in  the  stationary 
hinge  locaticms.  These  new  locations  are  determined  by 

ed  =  Li[L2-(L|-L?)l/2]/2jca2  =  [(H  +  Hd)/H]e  d2=  1^+  l^  =  lf  +  li  (B.50) 


where  Li  and  L2  are  the  new  lengths  of  the  chords  AB  and  B,  as  shown  in  Figure  B-1. 

MODIFIED  SUBMUNITION  RESPONSE 

Let  us  consider  the  previous  example  of  a  submunition  loaded  over  its  entire  length  when  the 
loaded  length  is  reduced  according  to 

H  1 
H  +  Hd“2 

Equation  (B.50)  and  geometrical  relationships  lead  to 


Chord  AB  angle 

00 

39^ 

Chord  BC  angle 

<Po 

31® 

Subtended  angle 

Vo 

62® 

Chord  AB  length 

Ll 

2.94  cm 

Chord  BC  length 

L2 

4.90  cm 

Chord  AB  height 

zi 

1.51  cm 

Chord  BC  height 

z2 

4.20  cm 

(Tcwnbined  heights 

D  =  d 

5.71cm 

Breadth  of  chords 

B 

2.52  cm 

Chord  AB  mass 

mi 

3.09  m 

(cm) 

Chord  BC  mass 

m2 

5.89  m 

(cm) 

Chord  AB  inertia 

U 

2.37  m 

(cm2) 

Chord  BC  inertia 

12 

14.78  m 

(cm2) 

Final  AB  angle 

e 

7t/2 

Final  BC  angle 

9 

36.9® 

Final  AB  height 

Zl 

0 

Final  BC  height 

z2 

3.92  m 
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Final  combined  height 

D  =  z2 

3.92  m 

Final  breadth 

B=Li 

2.94  m 

The  initial  values  of  the  coefficients  in  (B.18)  are 

dti)2/de  257.7  m 

(cm4) 

0)2 

-  337.2  m 

(cm4) 

P 

93.0 

cm3 

P 

M 

0 

-91.4 

cm 

and  the  final  values  are 

dG)2/de 

300.2  m 

(cm"*) 

0)2 

- 103.0  m 

(cm^) 

P 

89.9 

cm3 

P 

-  135.5 

cm3 

M 

-62.7 

cm 

The  coefficients  in  (B.20)  are 

Ao=1.13  Al=-1.79 

Bq  =  0.36  cm-1 

Bl=-0.11cm*l 

Co  =  0  Cl  =  0.83  cm-l 

Do  =  0.35  cm-3 

D1  =-0.26cm-3 

and  the  final  value  of  a  is  a  =  0.541  radius. 

The  pressure  that  triggers  the  mechanism  is 

Ps  =  (^  'jMo  =  67MPa 


The  mavimiini  fluid  pressuTB  is  p  =  126  MPa,  and  the  pressure  applied  throughout  the  motion 
required  to  generate  the  fluid  pressure  is 


Ps  =  93MPa 


For  the  choice  of  integral,  (B.32)  or  (B.33),  we  obtain 

93  MPa  <  P  <  210  MPa 
P>  210  MPa 


y2<0 

y2>0 
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and  the  ideal  impulse  required  to  produce  a  fluid  pressure  of  p  =  126  MPa  when  motion  ceases  is 

lo  =  149  ktaps  (1.49  MPa»s) 

Rgure  B-3  shows  the  PI  relationship  for  this  example  as  a  dashed  curve;  the  full  line  curve  is  the 
relationship  when  the  submunition  is  loaded  over  its  entire  length  (same  curve  as  that  in  Figure  B-2).  At 
high  pressure,  the  impulse  has  to  be  increased  by  about  20%  due  to  the  reduced  loaded  length.  The 
pressure  asymptotes  for  both  loading  cases  are  conparable  (Ps  =  94  MPa). 

PRESSURE-IMPULSE  LOADING  RELATED  TO  FRAGMENT  IMPACT 

CtHisider  an  idealized  fragment  ccmsisting  of  a  steel  rectangular  plate  traveling  at  a  velocity  Vf  altmg 
the  normal  to  the  plate  surface,  the  normal  being  perpendicular  to  the  axis  of  the  cylindrical  submunition. 
We  treat  the  case  in  which  the  length  of  the  plate  equals  the  height  of  the  cylinder  so  that  impact  occurs 
over  the  entire  height  Let  the  mass  per  unit  height  of  the  cylinder  and  plate  be  Me  Mf  req)ectively. 

By  assuming  that  the  impact  is  plastic,  we  have  the  velocity  of  the  cylinder  and  mass  as 


The  impulse  I  per  unit  height  is 

I=McV  =  Mf(Vf-V)  (B.52) 

The  impulse  I  of  our  PI  relationship  is  given  by 

I  =  2b  I  (B.53) 

where  2b  is  the  width  over  which  impact  occurs  (Figiffe  B-1).  Similarly,  if  the  contact  force  per  unit 
height  is  P,  we  have 


P  =  2b  P  (B.54) 

Let  the  impact  time  be  T  and  the  ring  deflection  of  the  plastic  hinge  mechanism  be  zi;  at  this  deflection  the 
hinges  are  collinear. 

We  assume  that  the  deceleration  of  the  fragment  from  velocity  Vf  to  velocity  V  is  constant.  This 
assumptitm  is  consistent  with  the  concept  of  an  average  contact  pressure  used  in  the  equation  of  motion. 
We  now  have 

zi=i(V  +  Vf)T  (B.55) 
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PRESSURE  (kbar) 


IMPULSE  (bar-s) 


Figure  B-3.  Pressure-impulse  relationships  for  a  fully  and  partially  loaded  canister. 


where  the  ccmtact  duration  may  be  expressed  as  T  =  I/P  =  I/P. 

Equations  (B.51)  through  (B.55)  can  be  solved  for  the  fragment  mass  and  velocity  in  the  forms 


2bl 


and 


In  (B.56)  and  (B.57), 


where  pf  is  the  fluid  density. 


■  2  [(Pzi/I)  -  (2bI/Mc)] 


Vf=(2Pzi/I)-(2bI/Mc) 


Me  =  p  27tah  +  pfTca^ 


(B.56) 


(B.57) 


As  an  example,  consider  the  submunition  example  presented  in  the  "Submunition  Response" 
section  and  let  the  fluid  be  water.  Then 

Me  =  114.5  g 
zi  =  1.02  cm 
2b  =  4.(X)cm 


From  Figure  B-2,  a  PI  pair  is 

P  =  300  MPa  I  =  142  ktaps  (1 .42  MPa«s) 

and  the  load  duration  is  T  =  49  ps.  Equations  (B.56)  and  (B.57)  give  Mf  =  18.4  g/cm  and  Vf  =  368  m/s. 
If  the  height  of  the  submunition  is  15  cm,  the  fragment  mass  is  276  g. 

An  improventient  is  effected  by  using  the  value  of  the  angle  a  associated  with  the  values  of  P  and  I. 
The  submunition  deflection  during  the  pulse  time  T  is  then  zi  -  zi,  where  zi  =  li  cos  0  with  0  =  (ic/2)  - 
a.  In  our  example,  we  have  zi  =  0.354  cm  to  give  the  deflection  during  contact  of  zi  -  zi  =  0.667  cm. 
Substituting  zi  -  zi  for  zi  in  (B.56)  and  (B.57)  results  in  Mf  =  31.1  g/cm  and  Vf  =  232  m/s.  Again,  for  a 
submunition  of  15  cm  height,  the  fragment  mass  is  467  g.  The  common  velocity  after  impact  is  V  =  50 
m/s. 
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APPENDIX  C 

TABLES  OF  TEST  RESULTS 
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Table  C-1.  Impact  tests  on  Chemical  Submunition  1. 


89  mm  (3.5  in.) 


Impactor 


89  mm  (3.5  in.) 


Impact 

Orientation 


^  89  mm  (3.5  in.) 

9mm  .470g 

(0.375  in.)  ^ 


Gun  7  A,  C,  F 

V  =326m/s 
MV  =605kg*m/s 
1/2  MV2  =0.986x105kgm2/s2 

Severe  fracture  of  upper  casing, 
fuze  flattened,  multiple  fragments 


Gun  a  A,C 

V  =234m/s 
MV  =435kg-m/s 
1/2  MV2  =  0.508  x105  kg  m2/s2 

Cracking  of  upper  casing,  some 
fragments,  fuze  crushed, 
circumferential  cracking  of  upper  weld 


Gun  9  A,C  Oun  1.1.  A 

V  =198m/s  V  =366m/s 

MV  =  366  kg*m/s  MV  =345  kg-m/s 

1/2  MV2  =  0.363  xIO^  kg  m2/s2  1/2  MV^  =0.632x10^kg  m^/s^ 

Cracking  of  upper  case,  no  fragments,  Severe  axial  fracture,  large 
fuze  caished,  circumferential  cracking  fragment  removed  from  side 
of  upper  weld 


Gun  10 

V  =175m/s 
MV  =  324  kgnn/s 
1/2  MV2  =  0.284  xIO^  kg  m2/s2 

Fuze  crushed,  very  slight  cracking 
of  upper  weld. 


Gun  12 

V  =358  m/s 
MV  =170kg*m/s 
1/2  MV2  =  0.305  xIO®  kg  m2/s2 

Slight  erosion  of  fuze,  no 
cracking,  impactor  fragmented 
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Table  C-1.  Impact  tests  on  Chemical  Submunition  1  (Continued). 


89  mm  (3.5  in.) 


Impactor 

89  mm  (3.5  in.) 


Impact 

Orientation 


„  89  mm  (3.5  in.) 

~940q  -470g 

(0.375  n.) 


Gon14  F,D 

V  =336  m/s 

MV  =318kg*m/s 

1/2  MV2  =  0.535  x105  kg  m2/s2 

Severe  dent  under  impactorf,  fuze 
separated  from  canister,  diaphragm 
blown  out,  slight  cracking  of  upper  weld, 
impactor  fractured 

Gun  15  D 

V  =273  m/s 

MV  =257kg*m/s 

1/2  MV2  =  0.351  x105  kg  m2/s2 

Severe  dent  under  impactor,  fuze  bent 
toward  impact  side,  diaphragm 
ruptured,  no  fractures  in  canister 

Gun  13 

V  =204m/s 
MV  =191  kg»m/s 
1/2  MV2  =0.195x10®kgm2/s2 

Large  dent  in  the  side  under  impactor, 
fuze  bent  toward  impact  side,  no 
fractures  in  canister 


Gun  17 

V  =197  m/s 
MV  =  257  kg*m/s 
1/2  MV2  =  0.183  x105  kg  m2/s2 

Upper  body  fractured,  large 
fragment  removed,  axial  cracking 


A,C  0un18 

V  =197m/s 
MV  =94  kgmVs 
1/2  MV2  =0.093x105kgm2/s2 
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Table  C-1.  Impact  tests  on  Chemical  Submunition  1  (Continued). 


Impactor 


89mrn(35irU 

38  mm 

(1.5in.)^|H||[|| 

89  mm  (3.5  in.) 
mm  |||||||||||M  ..Qdf)  n 

Qr«m  89  mm  (3.5  in.) 

(0.375  in.) 

\/ 

► 

Gun  19  C 

V  =296  m/s 

MV  =  103  kg*m/s 

1/2  MV2  =0.128x105|(gm2/s2 

Crushing  of  shoulder  at  impact  point, 
slight  cracking  of  top  weld 

V 

/ 

Gun  21  Q 

V  =343  m/s 

MV  =  142  kg*m/s 

1/2  MV2  =0244x105kgm2/s2 

Crushing  of  shoulder  at  impact  point, 
cracking  of  weld  at  impact  point,  some 
leakage  through  crack 

Gun  20  C 

V  =204  m/s 

MV  =102kg»m/s 

1/2  MV2  =  0.126  x105  kg  m2/s2 

Crushing  of  shoulder,  cracking  of  weld 
at  impact  point,  slight  leaking  at  crack 

1 

Gun  24 

V  =112  m/s 

No  water 

Little  damage  to  target,  impactor 
fractured  at  top  weld 

Gun  32  F 

V  =113m/s 

Water 

Little  damage  to  target,  fuze 
broken  off,  impactor  ruptured 

Table  C-2.  Impact  tests  on  Chemical  Submunition  2. 


Gun  28 

V  =334m/s 
MV  =  200  kgnn/s 
1/2  MV2  =0.334  x105  kg  m2/s2 

Nose  slightly  crushed,  impactor 
fragmented 


gun  27  C,D 

V  =200  m/s 

MV  =119kg^s 

1/2MV2  =0.119x10®  kg  m2/s2 

Upper  canister  fractured  but  top 
not  removed,  diaphragm  ruptured 

Gun  40  D 

V  =200  m/s 

MV  =120kg^s 

1/2MV2  =0.120x10®  kg  m2/s2 

Diaphragm  ruptured,  upper  canister 
intact.  Welded  construction 

Gun  31  C 

V  =299  m/s 

MV  =92  kg»m/s 

1/2MV2  =0.138x10®kgmV 

Top  broken  off  at  weld,  shock  tube 
broken  off 

Gun  45  D 

V  =197  m/s 

MV  =118kgtn/s 

1/2  MV2  =0.116x10®  kg  mV 

Diaphragm  ruptured, 

15%  ullage 

Gun  30 

V  =259  m/s 

MV  =78  kg«m/s 

1/2  MV2  =0.101x10®  kg  mV 

Side  dented,  no  fractures 

Gun  29 

V 

=  201  m/s 

MV 

=  61  kg*m/s 

1/2  MV2 

=  0.061  X 10®  kg  mV 

Side  lightly  dented,  no  fractures 
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Table  C-2.  Impact  tests  on  Chemical  Submunition  2  (Continued). 
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Table  C-3.  Impact  tests  on  scaled  Chemical  Submunition  2. 


I 


Fuze  removed, 
diaphragm  ruptured 


Gun  61  D 

V  =217  m/s 
MV  =15.8kg*m/s 
1/2  MV2  =0.0172x105  kg  m2/s2 

Diaphragm  ruptured 


Gun  62  D 

V  =205  m/s 
MV  =15.0kg»m/s 
1/2  MV2  =  0.0153  X 10®  kg  m2/s2 

Diaphragm  ruptured 


Gun  59 

V  =168  m/s 
MV  =12.3kgTn/s 
1/2  MV2  =  0.0103  X  10^  kg  rr?/$^ 


Gun  60 

V  =166  m/s 
MV  =  12.1  kgnrVs 
1/2  MV2  =  0.0101  X 10®  kg  n?l!? 


V  =198  m/s 
MV  =1.8kgmVs 
1/2  MV2  =  0.0018  X 10®  kg  m2/s2 

Diaphragm  ruptured, 
fuze  bent  forward 


V  =260  m/s 
MV  =1.3kg«m/s 

1/2  MV2  =0.0010x10®  kg  m2/s2 

Diaphragm  ruptured, 
fuze  bent  forward 


V  =197  m/s 
MV  =1,0kg*m/s 
1/2  MV2  =  0.0010  X 10®  kg  ni2/s2 

Fuze  bent  forward  slightly 


tttttttt 


Table  C-4.  Tests  to  define  P/I  curve  for  Chemical  Submunition  2. 


Orientation 


Impulsive  Loading 
Sheet  Explosive 


Quasi  -  Impulsive  Loading 
Body-to-Body  Impact 


Quasi  ■  Static  Loading 
MTS  Machine 


Uii 


Test  43 

I  =  91  ktaps 

C.D 

Upper  weld  fractured, 
diaphragm  ruptured 

Test  44 

I  =  63.7  ktaps 

F,D 

Fuze  removed, 
diaphragm  ruptured 

Test  42 

I  =  46.8  ktaps 

Cannister  slightly  damaged, 
diaphragm  and  fuze  intact 

I 

V  =230  m/s 
MV  =  90.9  kg*m/s 
1/2  MV^  =  0.104  X 10®  kg  m^/s^ 

Fuze  removed, 
diaphragm  ruptured 


V  =158  m/s 
MV  =  62.4  kg«m/s 
1  /2  MV^  s  0.049  X 1 0®  kg  m^/s^ 

Fuze  removed, 
diaphragm  intact 


Gun  55 

V  =155  m/s 
MV  =47.1kg-m/s 

1/2  MV2  =0.037x10®kgm2/s2 

Fuze  removed 
V^  =  76.6  m/s 


Gun  58 

V  =133m/s 
MV  =  53.3  kg»m/s 
1/2  MV2  =  0.035  x105  kg  m2/s2 

Fuze  removed 


Gun  50 

V  =119m/s 
MV  =47.1kg«m/s 

ieMV2  =  0.028  xIO®  kg  m2/s2 


Gun  54 

V  =111  m/s 
MV  =  32.9  kg«m/s 
1/2  MV2  =0.018x10®kgtn2/s2 

Slightly  dented  canisters 


Static  Test  3  -  Fuze  end  load 
Peak  load  =35,000  lb 
Peak  pressure  over  cannister 

=  28,520  psi  =  197  MPa 

Cannister  bulges  but  never 
ruptures 


Static  Test  1  -  Full  length 
Peak  load  =45,000  lb 
Peak  pressure  over  loaded 

area  =  19,000  psi  =  131  MPa 


Diaphragm  ruptured 


Static  Test  2  -  Unsupported  length 
Peak  load  =19,000  lb 
Peak  pressure  over  loaded 

area  =  13,500  psi  =  93.3  MPa 


Note:  All  tests  performed  at  half-scale. 
*1  tap  =  1  dyne  •  s/cm^ 


APPENDIX  D 

EXPLOSIVELY  ACCELERATED  PLATE  TESTS 


A  technique  was  devised  using  an  explosive  driver  to  accelerate  a  plate  to  impact  target 
submunitions.  Figure  D-1  shows  the  experimental  setup.  A  large,  dead,  soft,  aluminum  plate  (6061-T0) 
is  supported  on  top  of  a  steel  table  by  means  of  rigid  foam  collar.  The  steel  plate  is  lightly  attached  to  the 
bottom  of  the  aluminum  plate  and  lined  up  with  a  circular  hole  on  the  table  top  that  is  large  enough  to  pass 
the  steel  plate  through,  but  small  enough  to  stop  the  aluminum  plate.  Detonation  of  the  explosive  charge 
on  top  of  the  aluminum  plate  projects  the  steel  plate  at  the  submunition  located  under  the  table,  while  the 
dead  soft  aluminum  plate  forms  a  stretchable  diaphragm  that  covers  the  hole,  preventing  the  detonation 
products  from  trailing  the  steel  plate  and  obscuring  high-speed  photographs.  In  some  tests,  a  patented 
dilute  explosive  (developed  at  SRI)  that  has  an  adjustable  detonation  pressure  was  used.  The  use  of  the 
dilute  explosive  ensured  that  the  high  pressures  produced  by  conventional  explosives  would  not  break  up 
the  steel  plate.  Table  D-1  gives  the  tests  parameters  for  this  technique. 


D-1 


15.24  cm 


- 8.89  cm - 

-■ - 10.16  cm - - 

Figure  D-1 .  Experimental  setup  for  explosively  accelerating  steel  plates. 


Table  D- 1 .  Test  parameters  for  explosively  accelerated  steel  plates. 


Test  Number 

Thickness  (a) 
of  Aluminum 

_ (SIDi _ 

Thickness  (d)  Thickness  (h) 
of  Detasheet  of  DET 

(cm)  (cm) 

Thickness  (s) 
of  Steel  Disk 
(cm) 

Mass  of 
Steel  Disk 
(kg) 

Velocity  of 
Steel  Disk 
(m/s) 

TMD  Cal-1 

NONE 

0.11 

2.54 

0.98 

0.48' 

N/A 

TMD  Cal-2 

NONE 

0.11 

2.54 

0.98 

0.47 

N/A 

TMD  Cal-3 

0.64 

0.11 

2.54 

0.98 

0.48 

N/A 

TMD  Cal-4 

0.64 

0.11 

2.54 

0.98 

0.48 

N/A 

TMD  Cal-5 

1.27 

0.11 

2.54 

0.98 

0.48 

N/A 

TMD  Cal-6 

0.64 

0.11 

2.54 

0.98 

0.47 

N/A 

TMD  Cal-7 

0.64 

0.11 

2.54 

0.98 

0.48 

185.6 

TMD  Cal-8 

0.64 

0.11 

2.54 

0.98 

0.48 

180.4 

TMD  Cal-9 

0.64 

0.11 

2.54 

0.98 

0.48 

179.2 

TMD  Cal-1 0 

0.64 

0.11 

2.54 

0.98 

0.48 

196.4 

TMD  Cal-1 1 

1.31 

0.11 

5.08 

0.98 

0.48 

244.7 

TMD  Cal-1 2 

1.31 

1.28 

NONE 

0.98 

0.48 

318.2 

D-2 


APPENDIX  E 

AIR  GUN  PERFORMANCE  DATA 


The  relationship  between  gun  pressure,  sabot  and  impactor  mass,  and  impactor  velocity  is  depicted 
in  Figure  E-1,  where  sabot  velocity  is  plotted  versus  the  ratio  of  gun  pressure  to  sabot  mass  for  two 
different  types  of  sabots.  A  mathematical  fit  to  these  data  can  be  expressed  by  the  relationship 

V  =  C(P/M)1^  (E.1) 

where  V  is  sabot  velocity  in  meters  per  second,  M  is  sabot  mass  in  kilograms,  and  P  is  the  gun  pressure  in 
Pascals.  Values  of  C  that  are  used  to  fit  the  data  vary  between  0.28  and  0.34.  For  the  combinations  of 
mass  and  velocity  on  this  program,  the  value  of  C  was  found  to  be  about  0.30.  Using  this  value  for  C,  it 
was  possible  to  predict  impact  velocity  to  within  3%. 


V 


E-1 


PRESSURE  (GPa)  /  MASS  (kg) 


Figure  E-1 .  Gas  gun  velocity  calibration  curve. 


DISTRIBUTION  LIST 

DSWA-TR-96-19 


DEPARTMENT  OF  DEFENSE 

BALLISTIC  MISSILE  DEFENSE  ORGANIZATION 
ATTN:  CDR  HAGGERTY 
ATTN:  DANWHITENER 
ATTN:  DOUG  SCHAEFER 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN:  G  WEBER 

DEFENSE  SPECIAL  WEAPONS  AGENCY 
2CYATTN:  ISST 

3  CY  ATTN:  CDR  KENNETH  HUNTER 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
ATTN:  DTIC/OCP 

FIELD  COMMAND  DEFENSE  SPECIAL  WEAPONS  AGENCY 
ATTN:  DRBALADI 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  J  V  REPPA 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  ERIC  REECE 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


ARES  CORP 

ATTN:  CHUCK  MARTIN 


BATTELLE  MEMORIAL  INSTITUTE 
ATTN:  C  ALEXANDER 


JAYCOR 

ATTN:  CYRUS  P  KNOWLES 


KAMAN  SCIENCES  CORP 
ATTN:  VERN  SMITH 


DEPARTMENT  OF  THE  ARMY 


KAMAN  SCIENCES  CORP 

ATTN:  DENNIS  JONES 


U  S  ARMY  SPACE  STRATEGIC  DEFENSE  CMD 
ATTN:  LIBRARY 


KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 


USASSDC 

3CYATTN:  R  BECKER 


LOGICON  R  AND  D  ASSOCIATES 
ATTN:  DGAKENHEIMER 


DEPARTMENT  OF  THE  NAVY 

NAVAL  SURFACE  WARFARE  CENTER 
ATTN:  CHARLES  R  ELLINGTON 


LOGICON  R  AND  D  ASSOCIATES 
ATTN:  ETOTON 

MAXWELL  LABORATORIES  INC 
ATTN:  DRGGURTMAN 


DEPARTMENT  OF  THE  AIR  FORCE 

USAF 

ATTN:  RON  HUNT 

DEPARTMENT  OF  ENERGY 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  GPOMYKAL 


SRI  INTERNATIONAL 
aCYATTN:  A  L  FLORENCE 
2  CY  ATTN:  C  M  ROMANDER 
2  CY  ATTN:  J  COLTON 
2  CY  ATTN:  J  H  GIOVANOLA 
2CYATTN:  S  W  KIRKPATRICK 

W  J  SCHAFER  ASSOCIATES  INC 
ATTN:  J  CHERNAULT 


> 


Dist-1 


